Three years Swift— BAT Survey of AGN: reconciling theory and 

observations?* 



D. Burlon^, 
burlon@mpe . mpg . de 
M. Ajello^, 
maj ellotSisiac . Stanford . edu 
J. Greiner^, A. Comastri'^, A. Merloni^'^, and N. Gehrels^ 

ABSTRACT 

It is well accepted that unabsorbed as well as absorbed AGN are needed to 
explain the nature and the shape of the Cosmic X-ray background, even if the 
fraction of highly absorbed objects (dubbed Compton-thick sources) substan- 
tially still escapes detection. We derive and analyze the absorption distribution 
using a complete sample of AGN detected by Swift-BAT in the first three years 
of the survey. The fraction of Compton-thick AGN represents only 4.6% of the 
total AGN population detected by Swift-BAT . However, we show that once cor- 
rected for the bias against the detection of very absorbed sources the real intrinsic 
fraction of Compton-thick AGN is 20^g %. We proved for the first time (also in 
the BAT band) that the anti-correlation of the fraction of absorbed AGN and 
luminosity it tightly connected to the different behavior of the luminosity func- 
tions (XLFs) of absorbed and unabsorbed AGN. This points towards a difference 
between the two subsamples of objects with absorbed AGN being, on average. 
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intrinsically less luminous than unobscured ones. Moreover the XLFs show that 
the fraction of obscured AGN might also decrease at very low luminosity. This 
can be successfully interpreted in the framework of a disk cloud outflow scenario 
as the disappearance of the obscuring region below a critical luminosity. Our 
results are discussed in the framework of population synthesis models and the 
origin of the Cosmic X-ray Background. 

Subject headings: X-rays: general - Radiation mechanisms: non-thermal - X- 
rays: observations - Galaxies: active 



1. Introduction 



It is well known that absorbed active galactic nuclei (A GN) are needed to explain the 
shap e of the Cosmic X-ray background (CXB) spectrum (e.g. iGilli et al.ll2007l : iTreister et al. 
20091 ). A large f raction of them is irideed detected in the shallow and deep < lOkeV X-ray 



surveys (see e.g. iBrandt fc Hasingerll2005l ). No netheless, a lar ge fraction of Compton-thick 
(N/f > 1.5 X 10^^ H-atoms cm~^, for a review see IComastrill2004j ) AGN still escapes detection. 
Because of their large absorbing column density, these source s contribute a ~10-25% fraction 
of the CXB emission (IGilli et al.l 120071 : ITreister et al.ll2009l . respectively) at 30 keV, but at 
the same time ar e expected to be fairly numerous represe nting up to ~ 30 % of the entire 
AGN population (IRisaliti et al.lll999l : IWorsley et al.l 120051 ) . The advent of sensitive all-sky 
surveys above 15keV (e.g. Swift-BAT and INTEGRAL) opened the possibility to detect 
these objects. Indeed, it is above 10-15 keV that part of the primary continuum emission 
pierces through the veil of Compton-thick material, making it easier for these objects to be 
detected. Despite this fact, the early results from the Swift and INTEGRAL surveys showed 
that the fraction of Compton-thick AGN is merely a 5-10% of the total AGN population 
(see lAjellol l2009l . and references therein). Thus it might be that the fraction of Compton- 
thick AGN is intrinsically smaller than previously determined. Nevertheless, it should be 
taken into account that even above > 15keV, i nstruments are b i ased against the det ection 
of logN// > 24 sources. Indeed as shown in e.g. iGhisellini et al.l ( Il994l ). lAjellol ( 120091 ) . 50% 
of the source flux (between 15-55 keV) is lost if logNjy > 24.5. The fact that three of the 
closest AGN ever detected (i.e. NGC 1068, Circinus galaxy and NGC 4945) are absorbed 
by Nh > 10^^ atoms cm ^ suggests that indeed this is the likely explanation. 

In this work, we present and discuss the results of the most complete - until now - 
survey of AGN in the local Universe (z<0.1) using data from the Swift-BAT telescope. 
We performed a detailed spectral study of the three years sample of Seyfert-like objects by 
combining the hard X-ray information with the soft X-ray observations realized by different 
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missions in the recent past. To this aim we extracted Swift-BAT spectra and spectra in the 
0.3-10 keV band using archival XMM-Newton and Swift-XBT data. In a handful of cases 
we requested and obtained target-of-opportunity (ToO) observations with Swift for objects 
without previous coverage at soft X-rays. The use of X-ray data in the 0.3-195 keV band 
allows us to constrain robustly all the source parameters (including the absorbing column 
density) . 

This paper is organized as follows. In §2 we present the sample and discuss how the 
joint spectral analysis was performed. The general properties of the 15-195 keV continuum 
of AGN, the stacked analysis for the absorbed, unabsorbed, and Compton-thick spectra 
are presented in §3. §4 presents the observed Nh distribution. We then evaluate, for the 
BAT survey, the bias against the detection of the most obscured AGN and estimate -for 
the first time- the intrinsic absorption distribution. §5 presents the anti-correlation between 
the fraction of absorbed AGN (relative to the whole population) and luminosity, while the 
luminosity functions of AGN are derived in §6. The results of these analyses are discussed 
in §7 while §8 summarizes our findings. In this work we use a standard cosmology (Hq = 70, 
qo = 0, and Qa = 0.73). 



2. The Swift— BAT sample and data analysis 
2.1. The sample 



The Burst Aler t Telescope (BAT; iBarthelmy et al.l |2005[ ) onboard the Swift satellite 
( iGehrels et al.l 120041 ). represents a major improvement in sensitivity for imaging of the hard 
X-ray sky. BAT is a coded mask telescope with a wide field of view (FOV, 120° x 90° 
partially coded) aperture, sensitive in the 15-200 keV domain. Thanks to its wide FOV and 
its pointing strategy, BAT monitors continuously up to 80% of the sky every day achieving, 
after several years of the survey, deep exposure in the entire sky . Results of the BAT survey 
dMarkwardt et al.l bood : l Aiello et al.l l2008at iTueller et all boioh show that BAT reaches a 
sensitivity of ~1 mCrat0 in 1 Ms of exposure. Being the BAT survey not a flux-limited 
survey, but rather a significance-limited one, it is important to address how the survey 
flux limit changes over the sky area. This is often referred to as sky coverage, that is 
the distribution of the surveyed function of limiting flux. Its knowledge is very 

important when perforr ning populat i on stu dies as the ones described in the next sections. 
The reader is referred to lAjello et al.l (l2008al ) for how to derive the sky coverage as a function 



^ImCrab in the 15-55 keV band corresponds to 1.27x10 ^^erg cm ^ s 
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of the minimum detectable flux Fmm- This is defined as the sum of the area covered to fiuxes 
f < F • 



TV 



^ ( ^ Fmin ) ^ ^ ^ j ! fi ^ Fmin iX) 

i 

where N is the number of image pixels and Aj is th e area associated to e ach of them. A visual 



representation of the sky coverage is reported in lAjello et al.l (l2009bf ) which shows clearly 
the good sensitivity of BAT. The survey, in our analysis (15-55 keV), reaches a limiting 
sensitivity of ~0.6 mCrab (7.3 x 10~^^ erg cm~^ s~^ ). We ex post checked that cutting the 
sample at 50 % of the complete sky coverage (i.e. at 1.1 x 10~^^ erg cm~^ s~^) does not affect 
significantly the findings of this work, since the 35 objects below this limit do not populate 
a particular Nh range. 

The sample used in this work is the collection of non-blazar AGN detected by BAT 
during the first three years, mo re precisely between March 2005 and March 2008. This 



sample is part of the one used in lAjello et al.l ( l2009bl ) which comprises all sources detected 



by BAT at high (|b|>15°) Galactic latitude and with a signal-to- noise ratio (S/N) exceeding 
5 cr. All the 199 sources which are identified as non-blazar AGN (e.g. Seyferts) constitute 
the sample used in this work and are reported (along with their properties) in the table 
at the end of the paper. Note that the main sample, from which the sub-sample of AGN 
is derived, comprises 307 objects of which only 7 are as of today without identification. 
The incompleteness or the parent population is thus 2.3%. We note that in our sample 19 
objects are classified as 'Galaxies' or 'Liners'. We believe these are normal AGN (e.g. Seyfert 
galaxies) for which an accurate optical classification is not yet available in the literature. This 
is based on the fact that the average redshift, luminosity and absorbing column density are 
respectively 0.03, 4.5xl0^^erg s~^ and lO^^cm"^. These values are in good agreement with 
the ones derived from the rest of the sample giving confidence to our hypothesis that these 
objects are AGN. 



2.2. Extraction of Swift-BAT spectra 



For eac h source in our samp le we extracted a 15-195 keV spectrum following the method 
described in lAjello et al.l (l2008cl ). Here we recall the main steps: the details can be found in 
the aforementioned paper. For a given source, we extract a 12 channel spectrum from each 
observation where the source is in the field of view. These spectra are corrected for residual 
background contamination and for vignetting; the per-pointing spectra are then (weighted) 
averaged to produce the final source spectrum. Thus, the final spectrum represents the 
average source emission over the time-span considered here (three years). Moreover the 
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reader is referred to lAjello et al.l (l2009al ) for a discussion about the accuracy of the spectra 
extracted with this method. 



2.3. Extraction of the soft X-ray spectra 

The goal of the present work is to obtain a rehable estimate of the intrinsic absorbing 
column density for the BAT AGN. 102 objects (out of the 199 AGN in our sample) have a 
reliable estimate of the absorbing column density present in the literature. The l arge majority 



(86) of these naeasurements comes from earlier results of the Swift survey (iTueller et al. 



20081 : IWinter et al.ll2008l . l2009aU bl) while for the rest (16) we used single-source publications 
(see the table at the end of the paper for details). 

For all the objects without a Nh measurement we used the available follow-up observa- 
tions performed by two observatories (i.e. Swift-XBT and XMM-Newton in 83 and 12 cases, 
respectively). XRT was used preferentially, while XMM-Newton was used in a handful of 
cases (i.e. when the detection significance by Swift-XKT was too low to constrain the spec- 
tral parameters and/or the Nh). Only in 2 cases we could not find any XRT or XMM follow 
up (i.e. IH 2107-097) or the available soft X-ray pointing was not deep enough to extract a 
spectrum of the source (i.e. [VV2003c] J014214.0+011615). 

For the filtering and spectra extraction we used xrtproducts only on Photon Counting 
Level 2 event files (grades 0-12) and the standard ftools of Headas v6. 6. 3 software, and SAS 
v9. 0. for Swif t-XBT and XMM-Newton observations respectively. We used Xspec vl2A.0ad 
( lArnaudlll996l ) to perform, for each AGN, joint spectral fits between the 15-195 keV and the 
0.3-10 keV data. Normally the spectra of sources detected (in the 0.3-10 keV band) with 
sufficient S/N were re-binned as to have a minimum of 20 counts/bin. In a handful of 
cases spectra were re-binned as to have 10 counts/bin, and consequently Cash in place of 
statistics was adopted. 

During the spectral fitting stage we took into account , regardless of the spectral model 
used (i) the local Galactic absorption ( iKalberla et al.ll2005l ) and (ii) a normalization factor to 
account for the different inter-calibration of the two instruments and for a possible variation 
of the source between the observation epochs. Small differences in the computed value 
of Nh might be present when comparing observations taken at different times. The case of 
NGC 7582 is discusse d in ^4.2, but in general we warn that variation s of the column densities 
are expected (see e.g. iRisaliti et al.ll2002l |2009[ iBianchi et al.ll2009al and references therein). 
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3. General properties 

We discuss in tlie following the general properties of the sample, focussing on the hard 
X-ray continuum emission. The joint analysis is considered in this section only to the aim 
of splitting the parent population in un-absorbed and absorbed sources. Fig. [1] shows the 
luminosity-redshift plane for all the AGN in the BAT sample. The fc— corrected Lx lu- 
minosities (not corrected for absorption, but see §3.2 and §4.2) were computed according 
to: 



where Fx is the X-r ay flux in the 15-5 5 keV energy range as listed in the table at the 



end of the paper (see lAjello et al.ll2008bl . for details about the flux determination), and Vx 
is the photon spectral index. Throughout this work, absorbed sources are those with an 
absorbing column density (Nh) larger than (or equal to) 10^^ atoms cm~^. It is apparent 
from Fig. [1] that obscured AGN populate more densely the low-luminosity/low-redshift part 
of the graph with respect to the high-luminosity/high-redshift part. A Kolmogorov-Smirnov 
(KS) test between the redshifts of the two populations of AGN (absorbed and unabsorbed) 
shows that the probability that both classes are drawn from the same parent population is 
~ 1.4 X 10-3. 

The lines reported in Fig. [1] represent the current limiting flux of the BAT survey 
(~ 7.3 X 10"^^ erg cm"^ s"^, solid line) and a much brighter flux of 10"^*^ erg cm"^ s"^ 
(dashed line). In the shallower case it is apparent that sources detected are preferentially 
absorbed with logNn larger than 22. 



3.1. Analysis of the 15-195 keV continuum 

We performed a simple power-law fit to the BAT spectra in the 15-195 keV band to 
derive a measurement of the photon index. These values are reported in the last table 
with the corresponding 90% uncertainties. The distribution of photon indices of all the 
BAT AGN is compatible with a Gaussian distribution with a mean of 1.95±0.02 and a 
sigma of 0.27±0.02. We analyzed separately the distribution of photon indices of obscured 
and unobscured sources respectively. These are reported in Fig. [21 We find that the two 
distributions appear to be different with the one of obscured sources showing a mean of 
1.92±0.02 and a sigma of 0.25±0.02 while the one of unobscured AGN displays a mean of 
2.07±0.03 and a sigma of 0.27±0.03. This is also confirmed by the Kolmogorov-Smirnov 
test which yields a probability of 3 x 10^'^ that the two distributions are drawn from the 
same parent population. The CT sources are not the main drivers of the distribution of 
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Fig. 1. — Luminosity-redshift plane for unabsorbed sources (blue circles) and for absorbed 
sources (red squares). The solid line represents the limiting sensitivity of the Swift-BAT 
survey for a source with photon index equal to 1.9, while the dashed line represents the 
sensitivity of the survey at a brighter flux, i.e. 10~^° erg cm~^ s~^. 
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Fig. 2. — Photon indices distribution of absorbed AGN (logN// >22, red line) and unab- 
sorbed AGN (logNj:/ <22, blue line). The photon indices used are the ones derived in the 
BAT band only (i.e. 15-195 keV band) and reported in the table at the end of the paper. 
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Fig. 3. — Photon spectral index as obtained from fitting with a simple power-law the BAT 
spectrum alone, versus the Nh found from the combined spectra in the 0.3-195 keV energy 
range. Different color/symbol identify absorbed and unabsorbed sources. 
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absorbed sources. Indeed eliminating the 9 most absorbed sources from the sample does not 
in turn introduce an appreciable difference in the distribution. The KS test returns a null 
probability of ~ 1 x 10~^ in this case. Thus, there is an indication, albeit marginal, that 
absorbed sources display, on average, a harder power-law continuum than unobscured ones. 



According to iHopkins et al.l ( l2009l ) radiatively inefficient accretion flows (RIAFs) make 
the intrinsic X-ray spectrum of an AGN harder. This may cause to incorrectly classify, at 
energies below <10keV, an AGN as obscured if only simple estimators (e.g. hardness ratios) 
are used. This is not however the case for the present work, because: 1) BAT is able to 
sample the intrinsic power-law spectral index independently of the level of absorption, 2) 
all sources have sufficient signal-to-noise ratio to correctly derive the absorption level using 
XMM-Newton or Swift-XKT data in conjunction with BAT ones. Thus, while RIAFs might 
certainly affect the intrinsic shape of the 15-195 keV continua, we bel ieve that the differe nces 



observed here are ascribed to orientation effects (as shown already in lAjello et al.ll2008cl ). As 
it can be seen in Fig. |3]we show the scatter plot of photon indices versus the absorbing column 
density and indeed there is a weak indication of a correlation between the two parameters (the 
Spearman's rank is -0.27, and null hypothesis probability P ~ 1.2 x 10~^). The low, negative 
rank correlation coefficient and the P value indicate that a chance correlation can be excluded 
at more than the >3a level. However if the sources with unconstrained Nh are excluded 
from the sample the correlation is not statistically significant any longer (P = 1.3 x 10"'^). 
A mild correlation is anyway expected, because of the contribution of the higher reflection 
component of type 1 AGN, in the low energy channel of the BAT. The "softening" effect 
that is introduced when fitting with a simple power-law is further discussed in the following 
section. 

We note that three AGN (i.e. Mrk 766, IRAS 05480+5927, Mrk 739) show a very soft 
BAT spectrum (e.g. photon index >3). For Mrk 766 we analyzed XRT and BAT data 
jointly and found out that the intrinsic power-law seems to be ~2.0 and that indeed a large 
reflection component is required (the data show also the presence of a soft excess). The 
large reflection component is what makes (very likely) the BAT spectrum softer. The BAT 
spectrum of IRAS 05480+5927 is quite noisy and very soft. Nonetheless by a joint fit with 
XRT the photon index is constrained to be ~1.8. A cutoff-powerlaw is statistically (ftest 
probability > 3a) better but in turn requires the cutoff energy to be at 18 keV (in the 12-30 
keV range, 3a contours). As for Mrk 739, a joint fit with XRT gives a slope of 1.7. We 
note that again a cutoff-powerlaw is statistically preferred with a cutoff energy of ~ 50 keV. 
Our results are unaffected by the change of these three BAT-soft spectra to the flatter value 
reported for the joint analysis. 



Stacked spectra analysis - In order to investigate the global spectral properties of AGN 



- 11 - 



Table 1: Best-fit parameters for tlie stacked spectra. Errors are 90% CL and parameters 
witliout an error were kept fixed during tlie fitting stage. Tlie columns report the value of 
the photon index (F), normalization of the reflection component i?, the cut-off energy 
and the absorption. 



SAMPLE 


# Obj 


r 


R 


Ec 


Nn 












(keV) 


(1022 cm-2) 




All 


199 


^■'"-0.39 


< 4.50 


> 80 




5.2/8 


All 


199 


, on+0 08 

l-»U_g Qg 


1 00+°-''* 

^■'^^-0.36 


300 




5.2/9 


Absorbed^ 


96 


1 74+0.07 
J^- '^-0.07 




300 




3.3/9 


Unabsorbed^ 


92 


1 71+0.10 
J^- ' J^-0.06 


1 90+1.12 
1-^J_1,00 


300 




3.6/9 


Compton-thick 


9 


1.80 




88^11 




9.6/9 



^ For absorbed and unabsorbed AGN we have assumed that the inclination angle between the normal to the disk/torus 
and the line of sight is 60 and 30 degrees respectively. The class of absorbed AGN includes all AGN which are absorbed 
but are not Compton-thick (e.g. 22<LogNH <24). We did not include the 2 AGN for which Njj could not be calculated. 



we performed a stacking analysis of the AGN in the BAT sample. The stacked spectrum of 
several sources is produced performing the weighted average of all the spectra. The weight 
is chosen to be the inverse of the variance of a given bin and it is exactly the same procedure 
used to extract the spectra of each individual source. The same stacking technique has 
been already appl ied with success to both th e study of Seyfert galaxies and galaxy clusters 
detected by BAT (lAjello et al.ll2008bll2009a[ ). This stacking technique is appropriate for the 
stacking of background-s u btract ed spectra generated by coded masks telescopes. As already 
reported in Ajello et al. (j2009a ). this stacking technique allows to determine the average 
properties of a source population. The stacked spectrum of the 199 AGN is not compatible 
with a simple power law (x^/dof = 22.62/10). This is due to a substantial curvature of 
the spectrum around 30keV (see Fig. H]). We found that an acceptable fit to the data 
{x^/dof = 5.2/8) is achieved when using a PEXRAV model. In this case we find that the 
best fit parameter for the slope is l.TStoag (error are 90% CL). Given the small dynamic 
range of the BAT spectrum (15-195 keV) is impossible to disentangle uniquely the reflection 
component and the cut-off energy. Indeed, our best fit shows that, at 90% confidence, the 
normalization of the reflection component is consistent with zero while the cut-off energy 
is bound to be > 80keV. The parameters of this best fit are reporte d in Tab . [H I n order 
to avoid this degeneracy we fixed the cut-off energy at 300 keV (see iDadinal l2008l . which 
reported an average cut-off of 300 keV for BeppoSAX sources). The best fit parameters 
(reported also in Tab. [1]) for the photon index and the reflectior i are respectively l . SOln'n s 



and l.OOllJ:^^, which are in good agreement with the findings of iNandra fc Pounds! (119941 ). 
The normalization of the reflection component is compatible with the presence of a reflecting 
medium which covers an angle of 27r at the nuclear source. 



In addition we generated stacked spectra for unabsorbed (logNn < 22), absorbed 
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Fig. 4.— Stacked spectrum of the 199 Swift AGN. The dashed hne is the best-fit PEXRAV 
model. Note the substantial curvature of the spectrum around ~30-40 keV. 
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(22<logNH < 24) and Compton-thick (logNn > 24) AGN. A simple power-law fit to the spec- 
tra of unabsorbed and absorbed AGN yields that the best-fit photon indices are 2.13±0.06 
and 2.00±0.06 respectively. This is found to be in agreement with what derived from the 
analysis of the photon index of the two distributions: i.e. on average unabsorbed AGN 
have steeper spectra than absorbed ones. The indices derived from the stacking analysis are 
slightly steeper than the average ones derived from the photon index distribution because 
the stacked spectra show a significant curvature which makes the simple power-law fit not 
the most accurate one (e.g. reduced ^ 2.0). We thus decided to fit the stacked spectra 
with a PEXRAV model. Also in the stacked spectra of absorbed and unabsorbed AGN the 
cut-off energy and the reflection component cannot be determined uniquely. We thus flxed 
the cut-off energy to 300 keV. A fit to the stacked spectrum of absorbed sources yields a 
photon index of 1.741qq7, and a refiection component of R=0.55lo'35. 

The fit to the stacked spectrum of unobscured sources yields a photon index of 1.711q;o6 
and a refiection component of R=1.23]!^J oo- "^^^ uncertainties are large, however these results 
(which are reported in Tab. [1]) seem consistent with the unified mode l which predicts a la rger 
refiection component for unobscured sources (for a discussion see e.g. lAjello et al.ll2008ai and 
re ferences therein). In addition, our findings agree with the modeling of obscured sources 



m 



Gilli et al.l ( 120071 ). where the refiection efficiency for high inclination angles (expected for 



obscured AGN in the unified picture) is lower (0.88 rather than 1.3) than the one assumed 
for unobscured ones. 



3.1.1. The 15-195 keV Spectrum of Compton-thick AGN 

Finally we also investigate, for the first time, the average spectrum of Compton-thick 
AGN. Our first goal is to determine an empirical model which describes the 15-200 keV 
emission of Compton-thick AGN reasonably well and then later to interpret the features of 
the spectrum. Thus, we started fitting the stacked spectrum of the 9 Compton-thick AGN 
with a simple power-law model. The best-fit photon index is 2.04±0.09, but because of the 
spectral curvature this model represents a poor description of the data (x^/'^o/=42.1/10). 
The fit improves {x^ /do f =16.1/9) if we use an absorbed power-law model. In this case the 
best-fit photon index and absorbing column density are respectively 2.48lo'^8 ~ 
4.7l^;g X lO^'^cm"^. As a last step we tried to fit the stacked spectrum with an absorbed 
cut-off power-law model. We fixed the photon index of the power law to 1.8 to avoid de- 
generacy among the parameters. This model provides a good representation of the BAT 
data [x"^ /do f =9.7 / 9). The column density is consistent with being Compton-thick (Nh = 
289l}3^ X 10^^ atoms cm~^) and the cut-off energy is 82^'^g keV. The results of this fit are 
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Fig. 5.— Stacked spectra of obscured AGN (22 <logNH < 24), unabsorbed AGN (logN// < 
22) and Compton-thick AGN (logNw > 24 ) compared to the Cosmic X-ray background 
(CXB) as measured by lAjello et al.l ( l2008bl ). Note that both the CXB spectrum and the 
spectrum of Compton-thick AGN have been rescaled arbitrarily. 
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summarized in Tab. [TJ However, we caution the reader this model (zphabs in Xspec) takes 
into account only photoelectric absorption and it is used only as a functional form to show 
that the average continuum of Compton-thick AGN is indeed very curved. Fig. \5\ shows the 
average spectra of unabsorbed, absorbed and C ompton-thick AGN and compares it to the 
general shape of the Cosmic X-ray Background (jAjello et al.ll2008bl ). 



The peak of the stacked spectrum of Compton-thick AGN (at zf^O) is at almost twice 
the energy of the peak of the CXB (see Fig. testifying that if Compton-thick AGN are 
responsible for part of the emission at the peak of the CXB then the bulk of the population 
should be at z ^ 1. This seerns to be in agreement with the prediction of population synthesis 



models (e.g. iGilh et al.l 120071 : iTreister et al.ll2009r ). 



We also adopt a different strategy to check the results of the analysis of the stacked 
spectrum of the Compton-thick AGN. Instead of producing a stacked spectrum, we performed 
a simultaneous fit to the 9 spectra. In doing this we use the same baseline model for all 
the spectra (e.g. a power law or an absorbed power law) with a free normalization constant 
as to take into account the different source fluxes. As a first test, we tried as before a fit 
with a simple power-law model. This fit yields a photon index of 1.94±0.07 (in reasonable 
agreement with what found before) and a x^/'^o/=148.1/98. We then tried fitting with an 
absorbed power-law model. The improvement in the fit, with respect to the previous one, is 
very large (e.g. x^/dof =118.6/97) and the F-test confirms that the probability of obtaining 
such improvement by chance is ~ 4 X 10"^ The best-fit photon index is 2.29l|j;J^ and the 
column density is 2.4'^^^ x 10^^ cm~^. A fit with a cut-off absorbed power law model (with the 
photon index set to 1.8 as before) produces a slight improvement (e.g. = do f =113.1/97) 
yielding an absorbing column density of 1.8!l^o 3 ^ lO^^ cm~^ and cut-off energy of 128^^^ keV. 
Within the statistical uncertainties, the results of this new analysis appear to be consistent 
with the results of the analysis of the stacked spectrum. 

Our next step was then to fit the stacked spectrum of Compton-thick AGN with a 



more physical model. In this case we adopted an improved version of the model of lYaqoob 
( 19971) which fully treats rela tivistic Compton sca ttering, i.e. the MYTorus model by 



Murphy fc Yaqoob (2009) and Yaqoob et al. (2010]E This model provides tables for the 
attenuation of the continuum emission (transmitted through the torus) and the scattered 
component comp uted via Monte Carlo simulations (a similar model can be also found in 
Matt et al.lll999l ). In principle both the transmitted and the scattered components should 
be fitted to the spectrum to ensure self-consistency of the model. In practice, because of the 
many model parameters and the limited energy bandpass of BAT, a fit with both compo- 



^The model is available at www.niytorus.com 
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Fig. 6. — Confidence contours (1, 2, and 3 a) for the photon index and the absorbing column 
density of the MYTorus model fitted simultaneously to all the 9 Compton-thick AGN spectra. 
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nents was not successful (e.g. /dof > 3). However in this first exercise the normahzation 
of the scattered component was a factor > 10 larger then the transmitted one. We then tried 
to fit the two components separately to understand whether one component is dominating 
over the other one. The best fit using the transmitted components yields a /do f =20.6/ 9 
and is still thus a fairly poor fit. Instead we achieved a good fit {x"^ /do f =11.3/9) using the 
scattered component alone and an orientation of the torus (with respect our line of sight) 
of 60 degrees. The best-fit photon index is 2.17loj? while the absorbing column density is 
Njif = 3.7lfg X 10^^ cm~^. We also attempt a simultaneous fit with the MYTorus model to 
all the 9 Compton-thick AGN spectra, leaving the normalization of each spectrum to be a 
free parameter of the fit. Again it appears that the scattered component is dominating over 
the transmitted one. Indeed we achieve a good fit to the data using the scattered component 
alone (e.g. =111.1/97) and an inclination of the torus of ~60 degrees. The best fit photon 
index is 2. 08 ±0.10 and the absorbing column density cm ^. Fig. [H] shows the 

confidence contours of these two parameters. If we remove the most obscured AGN from the 
fit (i.e. NGC 1068) the index and become respectively 2.05±0.10 and 2.7~!:li x 10^^ cm'^ 
showing that our results are not driven by just one particular source. 

From the best fit using the MYTorus model (either to stacked spectrum or the simulta- 
neous fit) we derive that only ~30 % of the intrinsic nuclear fiux is observed in the 15-155 keV 
band. Finally, we note that the results presented in this section do not change if we remove 
the two brightest objects in the CT sample, nor if we remove the most absorbed source (e.g. 
NGC 1068, see above) from the sample. Nevertheless, given the paucity of CT AGN in our 
sample and the fact that they span one dex in absorbing column density, the results of this 
section must be taken with care as they might turn out not to be representative of the entire 
population of CT AGN. 



3.2. Luminosity distribution and spectral properties 

Luminosities of the Compton-thin AGN have been calculated through Eq. [21 therefore 
without taking absorption into consideration. We also tested whether the modeling of ab- 
sorption introduced a bias in the distributions, even in the hard band sampled by BAT. To 
this aim we fitted all the sources with logNn > 23.5 taking Compton scatter into consider- 
ation (i.e. we used cabs*zwabs*pow in Xspec) and compared the resulting de-absorbed fiux 
distribution with the one tabulated at the end of the paper. The cabs model has nonetheless 
some caveats that should be stated clearly: it assumes a constant Compton cross section 
equal to the Thomson one, so that it fails in describing the spectral hardening of the trans- 
mitted component due to the decay with energy of the Klein-Nishina cross section. Therefore 
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is typically used for spectra in the 2-10 keV energy band. Moreover, it does not take into 
account scattering from material out of the line-of-sight. The KS null probability (0.6) shows 
that indeed no appreciable difference is introduced by using a simple model for Compton-thin 
sources. As for the 9 Compton-thick objects, the fluxes have been de-absorbed as described 
in detail in §4.2. Fig. [7] shows the distribution of luminosity in the energy range 15-55 keV 
for the complete sample of AGN along with the distributions for the absorbed and unab- 
sorbed AGN. The median values of the two subclasses read a logarithmic value of 43.2 and 
43.8, respectively. A Kolmogorov-Smirnov test between these two populations reads a null 
hypothesis probability of 6.6 x 10"^, if we neglect those producing the peak at logNH=20 in 
Fig. [101 The distance between the two populations is exacerbated when taking into account 
all the sources with logNH=20, since the KS probability drops to ~ 10~^. Hence already 
by comparing the luminosity distributions there is evidence that the two populations are 
unlikely belonging to the same parent population, and this difference is independent from 
the modeling of absorption. This behavior is further discussed in detail in ^ where the 
X-ray luminosity functions of the two AGN classes are derived. 

Photon index vs. Luminosity - In general, a correlation between the photon index and 
the accretion rate (expressed as the ratio of the bolome tric luminosity and the Eddin gton 
luminosity) has been confirmed in numerous studies (see llshibashi fc Courvoisierll2010l and 
references therein). The general interpretation is that higher accretion rates lead to an 
increase in the photon density above the disk. This implies in general a more efficient cooling, 
and consequently steeper spectra. In this section we test whether the hard X-ray power law 
index is correlated with the luminosity in the BAT range, even if we are aware that the hard 
X-ray luminosity alone is an inaccurate proxy of the accretion rate. As it can be seen in 
Fig. [8] there is no indication that luminosity and photon index are correlated, confirmed by 
the Spearman's rank coefficient and null hypothesis probability respectively equal to 0.19 
and ~0.01. This might be produced by the fact that the obje cts in our sarnple hav e a broad 
distribution of black hole masses (and Eddington ratio^). Winter et al. ( 2009a ) recently 
showed the absence (in the local BAT-selected AGN sample) on average, of this correlation 
between the 2-10 keV de-absorbed luminosities and F, as well as between a proxy of the 
Eddington ratio and the spectral index. The correlation between the 2-10 keV photon 
index and the luminosity w as found not to be significant in many papers (iReeves fc Turner 
2000l : iBianchi et al.ll2009a[). and the significance is g enerally low also in the w orks where it 



is claimed as real (jPai et al.ll2004l : ISaez et al.ll2008l ). It is worth noting that IWinter et al. 



■^By definition XecIcI = Lsoi/LEdd and LecIcI = ^T^GcMirip/ gt, G and c being the gravitational constant 
and the speed of light, M and mp the mass of the black hole and of the proton, ct the Thompson cross 
section. 
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Fig. 7. — Luminosity distribution in the 15-55 keV energy range for the complete catalogue 
(grey) and for the absorbed sources (i.e. log Nh> 22, red). We also show (blue) the distri- 
bution for the sources with logNn < 22. The purple line shows the luminosity distribution 
of the absorbed sources without correction for absorption. 



(120081 ) reported a positive correlation between the photon index and the 2-10 keV flux of 
individual sources, i.e. at high fluxes t he sources tend to show steeper spectra. I nterestingly, 
Swift-B AT beamed AGN (Fig. 2 in IChiseUini et~aD boiol : lAiello et all l2009bh do show a 
rather remarkable relation between the 15-55 keV luminosity and the spectral index. This 
relation is even stronger when looking at samples selected in the Fermi /LAT energy range. 
Nonetheless, the reason of the absence of this correlation in the objects of our sample is 
beyond the aims of this work. 

Hard X-ray flux versus Absorbing column density - Fig. [9] shows the distribution of the 
sources in the flux-Nn plane. The absence of objects under ~ 7 x 10~^^ erg cm~^ s~^ (which 
is represented by the horizontal line) reflects the sensitivity of Swift-BAT in the energy range 
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Fig. 8. — Photon spectral index resulting from the simple power law fit as a function of 
the BAT luminosity in the 15-55 keV range. The unabsorbed sources are shown with blue 
circles, while the absorbed sources are shown in red squares. The 9 Compton-thick sources 
are highlighted with crosses, and their de-absorbed luminosity is connected to the absorbed 
one by horizontal lines. 
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we selected. Indeed we see no correlation between these two parameters (Spearman's rank 
coefficient 0.04, and null hypothesis probability P = 0.64). Note that the absence of sources 
at small Nh and high fluxes, i.e. the top left region of the plot, reflects the tendency of the 
more luminous sources (which are intrinsically less numerous) of being unabsorbed. Indeed 
we showed in Fig. [1] that at a higher flux the sources are typically absorbed by columns in 
excess of 10^^. This explains the handful bright objects, i.e. with Fx > 10~^° erg cm~^ s~^, 
in the range 22 <logNH< 23. We discuss in further detail the nine Compton-thick objects 
in M 



4. Absorption in the local Universe 

4.1. Observed Nh distribution 

We derived the photoelectric absorption by the combined fit of the BAT spectra (aver- 
aged over 3 yrs.) and the available follow up observations in the 0.3-10 keV energy range 
(e.g. XMM-Newton or XRT). As described in §2.3 we adopted the value reported in the 
literature for ~ 50% of the sample, and the specific reference is reported in the table at the 
end of the paper. The absorption distribution is reported in Fig. [101 For 33 AGN (out of 197) 
the absorbing column density we measured was found to be consistent with (or smaller than) 
the Galactic absorption in the direction of the source. When we could not constrain the ab- 
sorption or when its value was consistent with the Galactic one, we put the value logNH=20, 
which in turn produced the high peak in the distribution in the bin 20<logNH<20.5. Note 
that there are a handful of cases taken from the literature in which the column density is 
found to be lower than 10^*^ cm~^. 

When considering the whole AGN population we find that 53±4% (1 a statistical error) 
are absorbed by column densities > 10^^ cm~^. We find that the number of objects whose 
Nh is greater then 10^^ cm~^ is 9/197. Thus the fraction of highly absorbed sources, know 
as Compton-thick AGN is 4.6^^5% (1 a statistical error), all already known in the literature 
as extremely absorbed sources. These objects are highlighted in the table at the end of 
the paper. The fraction of Compton-thick sources at ~ 10~^^ erg cm~^ s~^ was predicted 
by populati on synthesis models t o be, at the typical fluxes sampled by BAT in the range 



7-15% (see IComastri et al.ll2009l : iTreister et al.ll2009l . and references therein). 



An updated compilation of the most recent survey results in the hard X-ray band (e.g 
above lOkeV) is reported in Tab. [2J It is clear that most, if not all, of these re sults indicate 



a lack of Compton-thick AGN if compared to the ~30 % fraction found by iRisaliti et al. 



fll999h 
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Fig. 9. — Scatter plot of the 15-55 keV fluxes of the sources vs. colum n density. Blue stars 
and purple triangles represent respectively the AGN already present in Ixueller et al.l ( 2008 ) 
and other references (see last table). Sources treated separately in this work are shown as 
red dots (follow up by Swift-XKT) and green squares (follow up by XMM-Newton). The 
horizontal dashed line represents the sensitivity limit of the 3-year survey (note that we 
draw a dotted line in the CT regime, where also BAT is biased, as discussed in § 14. 2p . 
Absorbed and un-absorbed fluxes (connected by vertical dash-dotted lines) are drawn for 
the 9 Compton-thick sources. 



-23- 




log(NH [cm ']) 



Fig. 10. — Observed column density distribution. The peak at logNH= 20 is produced by 
the sources for which absorption was neghgible (of the same order of the galactic one) . 
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Table 2: Observed fraction of absorbed and Compton-thick AGN, relative to the whole 
population for different X-ray surveys. 



Reference 


# Obj 


Completeness 


% Absorbed 


% C-thick 


Band 
(keV) 


Instrument 


Markwardt et al., 2005 


54 


95% 


~64% 


~10% 


15-200 


Swift-BAT 


Beckmann et al., 2006 


36" 


100% 


~ 64% 


- 10% 


20-40 


INTEGRAL 


Bassani et al., 2006 


56' 


77% 


~ 65% 


~ 14% 


20-100 


INTEGRAL 


Sazonov et al., 2007 


91 


93% 


~ 50% 


~ 10-15% 


17-60 


INTEGRAL 


Ajello et al., 2008a 


24 


100% 


~ 55% 


<20% 


14-170 


Swift-BKT 


Tueller et al., 2008 


103 


100% 


~ 50% 


~ 5% 


14-195 


Swift-BAT 


Paltani et al., 2008 


34"= 


100% 


~60% 


<24% 


20-60 


INTEGRAL 


Delia Ceca et al., 2008 


62 


97% 


~ 57% 





4.5-7.5 


XMM-Newton 


Malizia et al., 2009 


79'* 


100% 


~43% 


7% 


20-40 


INTEGRAL 


Beckmann et al., 2009 


135*= 


~97% 


~ 44% 


~4% 


18-60 


INTEGRAL 


This work 


197 


100% 


~ 53% 




15-195 


Swift-BAT 



^ The complete sample is 42 AGN, 36 of which are Seyfert galaxies. 

^ The complete sample is 62 AGN, 6 of which are Blazars and 14 are unindentified. 

Since the Paltani et al. sample may contain a fraction of spurious sources, we restricted their s ample to a limiting 
significance of 6 cr. Above this threshold all sources are identified (see Tab. 2 in lPaltani et al.ll2008l) . 

There are 88 objects reported to be at significance > 5.2<j. 79 of those are Seyfert galaxies, the remaining being Blazars. 

The complete sample c omprises 187 objects w ith > 3cr significance in the 18-60 keV energy band. According to the 
authors (see Sect 4.1 in lBeckmann et al.ll2009i ) there are 135 Seyfert galaxies with measured absorption. Only 7/187 
sources are listed generically as AGN without information on redshift. 

4.2. The BAT bias and the intrinsic Nh distribution 

BAT is the least biased X-ray instrument, particularly when comparing it to 2-10 keV 
telescopes, for the detection of obscured objects. Nonetheless, even in the > lOkeV band a 
relevant fraction of the source flux might be lost if the source is Compton-thick. In order 
to show this effect we computed the ratio between the observed and the intrinsic nuclear 
flux of an AGN for increasing column densities. In this exercise, we took both photoelectric 
absorption and Compton scattering into account using MYTorus. The nuclear, intrinsic, 
emission has been modeled as a power law with a photon index of 1.9. The results are shown 
in Fig. [11] From this plot it is apparent that the BAT survey is unbiased up to logNn ~24 
and then becomes biased against the detection of Compton-thick objects. Furthermore, the 
absorption bias affects much more severely the 2-10 keV band already for logNn >23 (see 
same Fig. [11]) . 

We also performed joint spectral fits of XMM-Newton and Swift-BAT data for all the 
Compton-thick objects in our sample to determine the fraction of the intrinsic flux which 
is seen in the 15-55 keV band. The details of this analysis will be reported in a future 
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Fig. 11. — Ratio of observed to intrinsic nuclear flux for an AGN, with a power-law spectrum 
with index of 1.9, as a function of the column density of the torus as seen edge-on. The 
red solid line shows the ratio for the 15-55 keV band for an AGN in the local Universe 
(e.g. the typical case for the BAT survey). The dashed line shows the ratio for the 2- 
lOkeV band for an AGN at redshift ~ 1, i.e. the typical red shift for AGN detected in the 
XMM-Newton/ Chandra surveys (see e.g. lAkylas et al.ll2006[ l. Note that the BAT survey 
is unbiased (e.g. ratio ~1) for logNn <23.5. The black data points show the ratio of 
observed to intrinsic flux estimated for the BAT Compton-thick AGN taking into account 
only photoelectric absorption. The gray data points show the ratio also when Compton 
scattering is fully taken into account. 



publication (Burlon et al. in preparation), however the results are summarized in Tab. [3] 
and in Fig. [H Only for SWIFT J0601.9-8636, NGC 5728, CGCG 420-015, and NGC 1068 
we used the values available in the literature. Indeed, XMM-Newton data are not available 
for SWIFT J0601 . 9-8 636 and NGC 5728, s o we used the values derived with Suzaku by 



Ueda et al.l (120071) andlComastri et al.l (120101 ) , respe ctively. NGC 7582 was largely discussed 

Bianchi et al.l (l2009bf l. as for the variations of a factor ~ 2 in 



PiconceUi et al.l ( 120071 ) and 
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the thicker absorber. We used the XMM-Newton observation taken in 2005, and therefore 
tagged this object as Compton-Thick. Variations of the thick absorber can in principle take 
place for other sources as well, altering the fraction of CT objects according to the observation 
used for the analy sis. The Compton-thick nature of CGCG 420-015 (aka IRAS 04507+0358) 
was discovered by lSevergnini et al.l ( l2010l ). It was known to be a highly absorbed Sy2 galaxy, 
but with the use of deep (100 ks) Suzaku observations, in addition to the ones available, the 
authors found Nh to be ~ 1.3 — 1.5 x 10^^ cm~^, with minor variations according to the 
model used to fit the combined spectra. NG C 1068 is a c ompl ex object which has been 
analyzed in detail in the past. According to iMatt et al.l ( 120041 ). the nuclear emission is 
completely suppressed and the source is s een only in reflec ted light. The column density 
is probably in excess of 10^^ atoms cm~^ (IMatt et al.l 120041) and t he reflection component 
is of the order of a few % of the nuclear flux (llwasawa et al.lll997l ). For these reasons the 
position on the plot of NGC 1068 should be considered tentative. Fig. [TT] also shows the 
difference in the observed-to-intrinsic flux ratio when also Compton scattering is taken into 
account (using the MYTorus model). Note however that modeling with MYTorus contains 
an implicit assumption, i.e. that the scattering material has a toroidal geometry with given 
parameters (half-opening angle of 60°, corresponding to a covering factor of AQ/An = 0.5). 
The line-of-sight of the observer is a flxed parameter of the flt. 

It is interesting to note that there is relatively good agreement between the model line 
and the observations of single Compton-thick objects detected in the Swift-BAT survey. In 
Fig. in] we showed the fluxes of the 9 Compton-thick AGN, before and after the correction 
for the missed flux, connected by dash-dotted lines. The horizontal line representing the 
limiting flux of the survey was turned into a dotted line to visually indicate the increasing 
bias in the CT regime. 



Table 3: Ratio of the observed flux with respect to the nuclear (unabsorbed) flux in the 
15-55 keV band for the 9 Compton-thick objects in the sample. Errors are quoted at 90 % 
CL. Four sources were fltted using XMM-Newton and Swift-BAT . For the remaining objects 
we used values from the literature. 



Source 


Fobs /Fnucl 


Error 


Nh 
(10^4 cm-2) 


Nh error 
(10^4 cm-2) 


NGC 3079 


0.15 


0.07 


5.40 


(-0.65,6.10) 


NGC 3393 


0.30 


0.10 


4.50 


(-0.36,0.62) 


SWIFT J0601. 9-8636 " 


0.94 


0.15 


1.01 


(-0.38,0.54) 


NGC7582 


0.97 


0.25 


1.10 


(-0.05,0.05) 


NGC 5728 >> 


0.82 


0.15 


1.0 


(-0.1,0.1) 


NGC 6240 


0.69 


0.25 


1.83 


(-0.23,0.22) 


NGC 1068 = 


0.02 


0.01 


>10 


(-5, >10) 


NGC 424 


0.92 


0.56 


1.99 


(-0.40,0.32) 


CGCG 420-015 


0.85 


~0.20 


1.46 


(-0.11,0.07) 



^ lUeda et al.l ||2007| ') '^ IComastri et al.l l l201(tl '^ iMatt et al.1 l l2004h and llwasawa et al.l lll997l ') ISevergnini et al.l l l2010h 
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Essentially the absorption bias limits the detection of Compton-thick objects only to 
those with bright (intrinsic) fluxes and in the very local Universe. Indeed three of the most 
famous Compton-thick objects (NGC 1068, NGC 4945 and the Circinus galaxy) are also 
among the closest known AGN. Thus the distribution reported in Fig. [TU] compares sources 
detected at different limiting intrinsic fluxes. It is possible to correct for this effect by taking 
into account the selection effect due to the large column density. The intrinsic absorption 
distribution can be expressed as: 

where the diV/dLogNn is in unit of sr~^ per logarithmic bin of Nh , 5* is the observed 
source flux, and dA^/dS'(NH) is the logA^-logS* of sources in a given logNn bin. The minimum 
ohserved?lViX (Smin) of integration should be set so that the limiting intrinsic flux is the same 
for all the bins. In this way the absorption distribution derived is representative of the density 
of sources at the same limiting intrinsic flux. 

The relationship between observed and intrinsic flux can be expressed as S°^^ = K{Ni{) 
where ^^'(A^h) is the ratio plotted in Fig. [TTl Thus Smin can be set to 10^^^ x K{Nii) to 
produce a uniform absorption distribution for sources with an intrinsic flux greater than 
10~^^erg cm~^ s~^. Eq. [3] can thus be rewritten as: 



dN 



A{Nh] 



a 



(4) 



dLogNu 1 

where A{Ni{) and a are the normalization and the index of the log N-log S in a given logNn 
bin. Here we assumed that the source count distribution can be approximated with a power- 
law function (e.g. dN/dS = AS'"). This assumption is well justifled by the fact that the 
source co unt distribution of t he entire BAT sample is well represented by a single power law 



[e.g. see lAjello et al.ll2009b[) and that the sourc e count distribution of AGN shows a break 
at much lower fluxes (e.g. ICappelluti et al.ll2007l ). In order to derive the intrinsic absorption 



distribution we used 0.25 dex logNn bins. For each bin, a source count distribution has 
been constructed and a power law has been fltted to it employing a Maximum Likelihood 
algorithm. This power law was then used to estimate the density of objects (per steradian) 
in a given logNn bin (i.e. Eq. |4]). All the power-law indices were found to be compatible 
with a Euclidean index (i.e. 2.5 for a differential distribution). The index a has been flxed 
to 2.5 for all those bins with less than 5 objects. We checked that flxing the index to 2.5 for 
all bins or allowing it to vary does not change the results presented here. For each logNn 
bin the error on the density of sources is derived from the Poisson error on the number of 
sources present in that bin to preserve the original counting statistics. 
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Fig. 12. — Intrinsic absorption distribution compared to the observed one (dashed hne). 
Errors bars are derived from counting statistics in each bin. 



The intrinsic absorption distribution is shown in Fig. [T2j From this we derive that the 
intrinsic fraction of Compton-thick sources {logNn >24) is 20+^ %, where the error is only 
statistical. It is apparent that the observed distribution starts to deviate from the observed 
one only for logN// >23.5. We also derive that the intrinsic fraction of absorbed sources is 
65±4 %. A systematic error in estimating the intrinsic fraction of Compton-thick sources 
might arise from both the power-law indices fitted in Eq. H] and from the K{Nii) correction 
factor (e.g. the curve plotted in Fig. [TTI) . Fixing all the a parameters to 2.5 or allowing 
all of them to vary changes the above fraction of about 1 %. Thus the exact shape of the 
logA^-logS* in each logNn bin does not contribute a large systematic error. 

On the other hand, the knowledge of the fraction of transmitted flux (e.g. K{N}i) and 
Fig. [TTI) plays a major role in the derivation of the density of Compton-thick AGN. The 
angle at which our line-of-sight intersects the torus and the power-law index of the intrinsic 
AGN spectrum can modify the fraction of Compton-thick AGN. Playing with these different 
parameters we derive that the systematic uncertainty on the fraction of Compton-thick AGN 
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is ~5 %. 

Also iMalizia et al.l (|2009[ ). using INTEGRAL, showed that the fraction of Compton- 
thick AGN is hkely larger than the observed ~5%. Instead of correcting for the missing 
population, they adopt a redshift cut (z<0.015) which would ensure, according to the au- 
thors, to have a complete sample. In their sample of 25 AGN, they found 6 Compton-thick 
AGN, thus the fraction of Compton-thick AGN is 24^g^ % of the total population in agree- 
ment with our estimate. 

As visible from Fig. [TTl there is also a slight overestimate (although compatible within 
1 o" with the observed density) of the intrinsic density of objects with LogN// ^ 20. This 
is due to the fact that a few objects in that bin have actually a lower column density that 
could not be effectively constrained in the 0.3-10 keV energy band. For all those objects 
the LogN// ~ 20 can be considered an upper limit to the true absorbing column density. 
Because of this, the source count distribution in the LogNj^ ^ 20 bin tends to overestimate 
the true intrinsic density. However, as seen in the Fig. [11] this effect is small. 



5. Anti-correlation of Absorption and Luminosity 



According to the AGN unified model ( lAntonuccil Il993l : lUrry fc Padovanil Il995[ ) all the 
different properties of AGN can be ascribed solely to orientation effect. Thus one should 
not observe variations of any other property with e.g. lumin osity, accretion rate and red- 
shift. However, already 30 years ago. [Lawrence fc ElvisI ( 1l982l ) reported the discovery of the 
anti-correlation of the fraction of obscured AGN (relative to the whole population) 

and luminosity. More recently different auth ors addressed the same issues with contra- 



dicting results. For example some studies (e.g. iTreister fc Urrvll2006l : iLa Franca et al.l 12005 



Delia Ceca et al.ll2008: Winter et al 



X-ray luminosity, while 



2009a 



Dwe 



Brusa et al. I I2OI0I ) suggest that F% decreases with 



Iv &: Pagg (2006 ) point to F % being independe i it of Ly. Also 
Sazonov fc RevnivtsevI ( l2004l ). ISazonov et al.l (120071 ) . and iBeckmann et al.l ( 120091 ) pointed 



out that an anti-correlation of the fraction of absorbed sources and X-ray luminosity seems 
to exist. 



Fig. [13] shows how the fraction of obscured AGN (those with logNn > 22) changes 
as a function of X-ray luminosity (in the 15-55 keV band) in our survey. The width of 
each bin has been chosen so that the number of object per bin is constant (~ 33). The 
errors on the number of absorbed sources and the total number of sou rces per bin have been 
propagated with the recipes for Binomial statistics (see lGehrelslll986l . in particular Tab. 6). 
Binomial statistics apply specifically when dealing with ratios of small numbers. The data 
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are correlated, as the Spearman's rank correlation coefficient is r^ = -0.94. The probability 
of a chance correlation (4.8 x 10~^) shows that the correlation is true at the 2.9a confidence 
level. 



We also show for comparison the INTEGRAL results of iBeckmann et al.l ( l2009l ). by 
converting the INTEGRAL 20-100 keV luminosities into 15-55 keV luminosities and selecting 
only non-blazar sources with z<0.3. The two data-sets are in reasonable agreement although 
the lower fraction of absorbed AGN in the INTEGRAL sample (with respect to the BAT 
sample) might be ascribed to completeness issues of the former. 

In the BAT sample, the absorbed AGN fraction displays a fiattening in the low luminos- 
ity regime, and the correlation becomes clearly non-linear {x^/dof > 3.5). Thus, we tried to 
fit the fraction of absorbed AGN with the empirical function proposed by iGilli et al.l (120071 ) 
of the form: 



F%{L 



X, 



{-Lx/Lc) 



Rhigh\^ 



(5) 



where Riow is the low luminosity asymptotic behavior, Rhigh the high luminosity one, and Lc 
is a "critical" luminosity at which the drop occurs. We fitted this function to the data: the 
best fit values for Riow, Rhigh, and log(Lc) are respectively 0.8, 0.2, and 43.7. This fit yields 
a x^/dof = 1.3 which shows that this fit is better than a simple linear relation between the 
obscured AGN fraction and (the logarithm of) luminosity. This result will be discussed in 
details in in § [61 



In the literature (e.g. iHasingerl l2008l : iDella Ceca et al.l 120081 ) the minimum luminosity 



considered for an AGN in the 2-10 keV energy range, was typically greater than 10 erg 
s~^. Since the low luminosity tail of our distribution (see Fig. [7]) extends to ~ 2.8 x 10^*^ 
erg s~^ we checked, for the sake of completeness, if by introducing a cut at 10^^ erg s~^, 
changed our findings. When transforming the 2-10 keV luminosity of 10^^ erg s~^ to the 
15-55 keV band, we find that only 7 sources fall below this limiting luminosity. The level of 
correlation of the data remains unchanged, while instead the x^/dof value for a linear fit, 
decreases to ~ 1.8. These results are shown in Fig. [13] (b). It is still apparent that the hnear 
fit (which shows a slope of —0.26 ± 0.05) is not a good representation of the data since the 
first bin and the last two show a fiattening of the fraction of absorbed AGN. Finally, it is 
worth consider ing that (i) the con t ribution from the ste llar population hardly extends above 
~ 10^^ erg s~^ (iRanalli et al.ll2003l : [Norman et al.l 120041 ). and (ii) in the 15-55 keV band this 
contribu tion is expected to d ecrease to even a lower fraction of the "bolometric" luminosity 
(see e.g. 



Voss fc Aiellol[2010[ ). 
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(a) The red dashed Une represent the fit with Eq. \5\ 
Grey stars rep r esent , for comparison, data from 



Beckmann et alJ (|2009[ ). Luminosities measured by IN- 



TEGRAL have been converted into the 15-55 keV en- 
ergy range. 
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(b) Same as above (a) but with a cut in luminosity, as dis- 
cussed in the text. The dashed (dotted) line represents 
a linear fit to the data (Icr uncertainties). 

Fig. 13. — Fraction of the obscured (i.e. Nh > 10^^ cm~^) AGN as a function of 15-55 keV 
luminosity (red squares). The data have been grouped such as to have an equal number of 
sources per luminosity bin. Errors on fraction of absorbed AGN have been computed using 
Binomial statistics (see lGehrelslll986[ ) and are drawn at la level. 
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Fig. 14. — Luminosity function of absorbed (logNjif >22, squares) and unabsorbed AGN 
(logNji/ <22, circles) with the respective best-fit models (solid and dashed lines). The lumi- 
nosity functions have been multiplied by to highlight the different positions of the break 
luminosities. 
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6. Luminosity Functions of AGN 



We estimate the X-ray luminosity function (XLF) of AGN using tlie 1/Vmax method 
(equivalent in our formalism to the 1/Va method). In this framework the XLF can be 
expressed as: 

^ dN 1 A 1 



dL 



X 



A Lx ^ Vmax,i 



where Vmax,! is the maximum comoving volume associated with the ith source. This is 
defined as: 



V 



z=zmax 



2=0 



dV 

dz 



Q{Li, z)dz, 



(7) 



where dV/d z is the comoving volume element per unit redshift and unit solid angle (see e.g. 



iHoggI 119991 ) ■ zmax is the maximum redshift after which the sources would not be detected 
anymore in the BAT survey and Q{Li, z) is the sky coverage of the survey. The XLF of the 
two different classes of AGN, obscured and unobscured, is reported in Fig. [1] 



We made a Maximu m Likelihood fit to the two different da tasets using a broken power- 
law of the form (see e.g. lUeda et al.ll2003l : iHasinger et al.ll2005l ): 



dN 



A 



dLx ln{10)Lx 



Lx 



71 



72' 



(8) 



The ML estimator can be expressed as: 

(^{Lx^i,Zi)V{Lx^i,Zi 



-2 5: In 



j^{Lx^z)V{Lx^z)dLx 



(9) 



Table 4: Best-fit parameters of X-ray Luminosity Functions in the 15-55 keV band 



SAMPLE 


# Objects 


Norm.-'- 


L-^ 


71 


72 


ALL 

ABSORBED 

ABSORBED^ 

UNABSORBED 


199 
105 
99 
92 


1.53e-5 
2.59e-5 
3.93e-5 
1.90e-6 


o.53tj;;i^ 

2fi+°"* 

2fi+° " 

1 Q4+0.48 


0.58i«:i^ 

o.5ilg:S 


Z.OU_Q 20 
+0.34 
^- ''^-0.30 
o ^■0+0.38 
^■DJ_0.31 

^■°°-0.31 



^ Normalization of the XLF expressed in units of erg~^ sMpc~^. 
In units of 10** erg s~^. 

^ Intrinsic XLF of absorbed AGN. The luminosity of the absorbed AGN with LogN^f > 23.5 have been de-absorbed with 
the method described in § 14.21 

The best-fit parameters are obtained by minimizing C Their 1 a error are computed by 
varying the parameter of interest, while the others are allowed to fioat, until an increment of 
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A£=l is achieved. This gives an estimate of the 68 % confidence region for the parameter of 



interest (lAvnilll976l ). The hkehhood function does not depend on the normahzation A since 
it cancels out in Eq. O Once the slope a is determined, the normalization is derived as the 
value which reproduces the number of observed sources. An estimate of its statistical error 
is given by the Poisson error on the number of sources used to build the XLF. 

The results of the ML fits to the XLFs of whole population of AGN and obscured and 
unobscured subclasses are summarized in Tab. IH We focussed mainly on the difference 
between the two subsamples of absorbed and unabsorbed AGN, but we used the total XLF 
in Fig. [16] in order to account for the obscuration- luminosity relation. It is apparent that 
the XLFs of the two classes of objects are not the same. In particular the 'break' luminosity 
L* is different at ^2.8( 7 level, with absorbed A GN having on average lower luminosity than 



unabsorbed ones. Also iDella Ceca et al.l ( 120081 ). analyzing a small sample of XMM-Newton 



AGN, found different XLFs for obscured and unobscured sources. However, in their case 
they cannot allow (presumably due to the low number of sources) L* to be a free parameter 
of the fit. In our case, this can be done and there is evidence (albeit marginal) that the 
typical luminosity of absorbed and unabsorbed AGN is different. The difference between the 
luminosity functions of absorbed and unabsorbed objects is however not a surprise. Indeed 
it is expected in view of the anti-correlation of the absorption fraction and luminosity (e.g. 
Fig. [T3|) . The two luminosity function are equal at a luminosity of 4x10^^ erg s~^. This is 
exactly the luminosity at which the fraction of absorbed objects is 0.5 (see Fig. [T3|) . For 
the very first time this trend is clearly seen in the luminosity function of absorbed and 
unabsorbed objects (as derived from the same energy band). 

We tested if this 'trend' of the break luminosity holds for smaller subsamples, namely 
if the very most absorbed sources show a significantly smaller value of L* . This exercise 
needs of course a sufficient number of sources per subsample, in order to constrain the 
parameters of the fit with the broken power law. Therefore we divided the parent population 
into four bins of absorption containing ~ 50 sources each. For each of the subsamples we 
computed the best fit values of the XLF as described before. In Fig. [15] we showed the break 
luminosity (in units of 10^^ erg cm^ s~^) versus Nh. This exercise confirms, albeit with a 
statistical significant uncertainty, that absorbed AGN appear to be intrinsically less luminous 
than unabsorbed AGN. We also performed an additional test in order to exclude that this 
finding is partially driven by the bias against the detection of the most absorbed AGN. We 
computed the intrinsic XLF of absorbed AGN by de-absorbing the AGN luminosities using 
the model described in §4.21 As clearly seen from Fig. [T^l this correction is negligible for all 
AGN with LogNH<23.5, modest for all those with 23.5 <LogNH<24 and relevant for AGN 
with LogNH>24. It has to be noted that the intrinsic XLF of absorbed AGN suffers from 
incompleteness at the lowest luminosities. Indeed, because of the effect of large absorption. 
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22 23 
LogN 



Fig. 15. — Break luminosity (L*, as in Eq. [8]) in units of 10^^ erg cm^ s~^, versus column 
density. The sample was divided into four equally populated subsamples of ~50 sources. 
Errors are quoted at 1 a level. 
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sources with an intrinsic luminosity large enough to be detected by BAT might be pushed 
below the BAT sensitivity. In order to avoid this problem we cut the sample at the minimum 
de-absorbed luminosity for which the K{Nh) correction (discussed in §4.2p was less than 1. 
This minimum, de-absorbed, luminosity is 2x10'^^ erg s~^ and we consider the BAT sample 
to be complete above it. The parameters of the intrinsic XLF of absorbed AGN are reported 
in Tab. HI It is clear that the intrinsic XLF of absorbed AGN is found to be in very good 
agreement with the XLF of absorbed AGN. This is because the bias against the detection 
of absorbed AGN is relevant only for LogNH> 24 and above this threshold the BAT sample 
contains very few objects. 

Finally, we test whether the anti-correlation of the fraction of obscured AGN and lu- 
minosity (found in the previous section) can be reproduced by the two different XLFs for 
absorbed and unabsorbed AGN. To this extent we divide the XLF of absorbed AGN by 
the XLF of the entire sample. We computed the 1 a error on this function via bootstrap 
with replacement employing 1000 bootstrapped samples. In order not to suffer from biases 
derived from the detection of Compton-thick AGN and also to avoid the need of correcting 
for them, we also excluded all objects with LogiVjy >24. As shown in the previous sections 
the BAT sample can be considered complete for all AGN with LogNn <24. The results of 
our analysis are shown in Fig. [16] along with the observed anti-correlation of the obscured 
AGN fraction and luminosity. The agreement within errors is clear at all luminosities and 
the decreasing trend in the absorbed fraction at low luminosity does not depend on the 
presence or the absence of Compton-thick AGN. The trend of the fraction of obscured AGN 
and luminosity can be perfectly reproduced by the analysis of the XLFs of the absorbed and 
unabsorbed AGN. In turn, this difference can be ascribed to the fact that on average ob- 
scured AGN appear to be less luminous. Even more interestingly, the ratio of XLFs predicts 
(in agreement with the observations) that the fraction of obscured sources decreases at low 
luminosities (i.e. Lx < 10"^^ erg s~^). However, we caution the reader that, due to the lack 
of sources, the obscured AGN fraction is, at low-luminosities, compatible within 3 a with a 
constant fraction of ~ 80 %. 



7. Discussion 



7.1. Compton-thick AGN and the Cosmic X-ray Background 



It is well understood that the shape and the intensity of the CXB cannot be reproduced 
successfully if Compton-thick AGN (logNn > 24) a re not invoked. They are required t o 
produce ~10% of the CXB intensity at ~ 30 keV JciUi et al.l[2007l : iTreister et al.lboOQl ). 
However, since the extreme absorption makes these objects faint at X-rays, Compton-thick 
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Fig. 16. — Fraction of the absorbed AGN versus luminosity in the 15-55 keV energy range. 
The hues are obtained by dividing the XLF of obscured AGN by the XLF of the whole 
population. The dashed (solid) line represents the whole sample (with a cut at 2 x 10*^^ erg 
s~^). The shaded bands represent the la uncertainty computed with a bootstrap analysis, 
blue for the whole sample and yellow for the sample with the cut in luminosity. 
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AGN have to be numerous (e.g. ~30 % of the whole AGN population) in order to fulfill the 
above requirement. Despite this general belief, all searches of Compton-thick AGN have so 
far highlighted a lack of Compton-thick AGN (see Tab. This work, which uses the largest 
AGN sample collected so far above 15 keV, shows that the detected {Taction of Compton-thick 
AGN is only 4:.6~^'l'}^%. At the fiuxes currently sampled by Swift -B AT popu l ation synthesis 
mo dels predict th i s fract ion to be either ~16% if one adopts the iGilli et al.l (120071 ) or ~8% 
the iTreister et al.l (120091 ) model. Thus our measurement appears to be subs tantially lower 
than both predictions, but compatible within ~2 a with what reported by iTreister et al. 
J2009h . 



However, one must take into account that even X-rays with energies larger than 15keV 
are absorbed if the source is mildly Compton-thick. Fig. [TT] shows how dramatic this effect 
is. Indeed, if the source has logNn ~24 then only ~50 % of the intrinsic, nuclear, flux is 
visible above 15keV. This fraction becomes much lower (e.g. a few %) if the source has 
logNn ~25. Thus the extreme low number of Compton-thick AGN detected in shallow 
surveys at hard X-rays is due to fact that only the population of intrinsically bright objects 
is being sampled. As such these objects are rare. For the flrst time, we use the knowledge 
of how X-rays are absorbed through a Compton-thick medium to determine the intrinsic 
absorption distribution. We flnd that Compton-thick AGN are 20lg% (statistical error) of 
the whole AGN population. Our result shows beyond doubt that Compton-thick sources are 
indeed intrinsically numerous, but due to the large absorption, very difficult to detect. 

The average spectrum of all Compton-thick AGN detected by BAT is, in the 15-195 keV 
band, compatible with a very curved spectrur n peaking around 4 0-50 keV (see Fig. H]). Since 
the peak of the CXB is located at ~ 25keV (lAjello et al.ll2008bl ). this implies that the bulk 
of the Compton-thick AGN must be located at zf» 1 in order to explain the missing fraction 
of the CXB. We employed the MYTorus model, that fully treats photoelectric absorption 
and relativistic Compton-scattering in a toroidal geometry, to fit the BAT spectra of the 
Compton-thick AGN. We derived that in the BAT band the spectrum of Compton-thick 
AGN is very likely dominated by the scattered component. The photon index is compatible 
with 2.0 while the absorbing column density is constrained to be Compton-thick using BAT 
data alone. We also found out that on average only ~30 % of the intrinsic nuclear fiux is 
visible in the BAT band. A more detailed an analysis on each of the Compton-thick AGN 
detected by BAT will be reported in a future publication. 

Summarizing we find that Compton-thick AGN: 1) are intrinsically as numerous as 
required by the AGN unified model , and 2) have a spectral shape which 'fits' the one required 
to explain the sh ape and the peak of the CXB. How ever, the predictions of population 
synthesis models (IGilli et al.l 120071 : ITreister et al.ll2009l ) in terms of observed Compton-thick 
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AGN are between a factor of 2 to 3 larger than what is currently observed by BAT. 



7.2. A receding torus or a Clumpy torus model ? 



It is generally accepted that all kinds of AGN are different manifestations of the same 
engine (i.e. the accreting super- massive black hole), seen under different line of sights. Thus, 
in the unified picture, obscuration depends solely on orientation effects. This hypothesis 
breaks down when showing that the fraction of obscured sources (relative to the whole 
population) decreases with luminosity. This result has been reported several times in the 
past, but it has been confirmed now by BAT (with a sample a factor ~2 greater than the 
one presented in lTueller et al.l (120081 )). 



Our results, like the previous claims of this anti-correlation, are inconsistent with the 
simplest predict ion of the unified model, for which F% should not depend on Lx- Already 
Lawrence! ( 1l99l[ ) tried to overcome this contradiction by proposing a receding torus model. 
In this model, the inner radius is set by the dust sublimation radius which increases with 
source luminosity. This model predicts that the fraction of obscured sources scales with 
luminosity as F% oc L~°-^. 



In basic terms (see iNenkova et al.l 120081 . for a detailed discussion), the inner radius is 
set from the following: 



Rd ~ 0.4 



L 



°' /1500 Ky-'p^ 



sub 



(10) 



where L is the bolometric luminosity and Tsub is the sublimation temperature of the dust. 
Given the angular dependence of F% (i.e. the angle at which a smooth torus becomes 
transparent to radiation), a constant H/R implies F% oc While attractive, this simple 

idea fails to reproduce correctly the dependence of the obscured fraction with luminosity, in 
particular when looking just at the hard X-ray portion of the spectral energy distribution. 



Indeed, in this work we showed that F% oc L 



-0.26±0.05 



for L 



X 



>10^2g^g g-l_ 



Nonethele ss, one has to c onsider that (i) the torus is likely clumpy, (ii) H/R is a function 
of luminosity ( ISimpsonI |2005| ). and (iii) that even if the dusty components of the torus do 
absorb optical, UV, and X-ray photons. X-ray obscuration can take place also in the -dust 
fr ee- inner region, in the pro ximity (< O.lp c) of the AGN. This has been largely discu ssed 



m 



Risaliti et al.l (120021120071 1. Recently fsee lHonig k Beckertll2007l : iNenkova et al.ll2008r i Fc 



>0 



was interpreted in the framework of a radiation-limited clumpy dust torus in which F% oc 
^-0.25^ which is much closer to our findings. In this scenario, the obscuration-luminosity 
relation is recovered in terms of probability that the photons intercept a sufficient number of 
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clouds in the line of sight. In addition, the work of iHasingerl (j2008[ ) s hows similar results in 
the 2-10 keV energy band, even if just proxies of Nh are used. Finally, iMaiolino et al.l (|2007| ) 
interpreted those early results by comparing them to a similar relation between the fraction 
of obscured AGN and the 0[III] luminosity (as well as other bands, but see their detailed 
description). According to the authors this was a signature of the variation of the covering 
factor of the AGN dust with luminosity. The variation of the coverin g factor of the torus is 



in ge neral- also associated with the "Iwasawa-Taniguchi" (IT) effect (llwasawa fc Taniguchi 
993h . i.e. the anti-correlation of the EW of FeKa an d X-ray l umino sity. Bianchi et al 



(l2009al ) discussed the agreement among the findings of iDadinal (120081 ) in the 20-100 keV 



regime, and their updated XMM-Newton sample. If the anti-correlation of F% and luminosity 
is interpreted as the variation of the covering factor of the toru s , then this can be in turn 
compared to the IT effect. From the results of iDella Ceca et al.l (120081 ). they found a slope 
of ^ 0.22, in v e ry go od agreement with the IT effect, with the slope (~ 0.18) presented in 
Maiolino et al.l (120071 ). and eventually with our findings (~ 0.26). Although these similarities 
are intriguing, they certainly deserve a profound investigation (which is beyond the goals of 
this paper). This is even more compelling when considering that the methods differ and in 
principle sample different materials. 

In the low luminosity regime (i.e. below 10^^ erg s^^), the behavior of the torus is 
likely more complex than what reported above. Different authors (e.g. 



Elitzur fc Shlosman 



20061 : lElitzur fc Hoi |2009|) showed that below a bolometric luminosity of ~ 10 erg s the 



torus obscuration region disappears. This results under the assumption that the clouds are 
gener ated by a disk- wind outflow, rath er then accreted fron a the g alaxy (iKrolik fc Begelman 
1988). In the framework proposed by Elitzur fc Shlosman (2006) and Elitzur fc Hoi ( 20091 ) 
the cloud mass outflow Md is proportional to the mass accretion rate, which is related to 
the bolometric luminosity via Mace = 0.02L4^/r]MQ yr~^, t] being the accretion efficiency. 
The ratio Md/Macc increases at progressively lower luminosities but since it cannot exceed 
unity, there is a limiting luminosity at which the system cannot sustain the cloud outflow 
any longer. 

At low luminosities we then expect absorption of X-rays to be less effective, due to 
the lack of obscuration in the torus. We should then observe a decrease in the fraction of 
obscured AGN at low luminosity. For the first time, and thanks to BAT, we are able to 
inspect the behavior of F% at low luminosities. Indeed, Fig. [16] show^ that the fraction 
of obscured AGN decreases at low luminosities as one would expect if the torus obscuring 
region would cease to exist. There are no other evidences (perhaps beside this one, but 



*We caution the reader that within 3a the behavior of F% at low luminosity is still compatible with a 
constant value of 0.8. 
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see Ivan der Wolk et al.l I2OIOI ) which shows that the torus disappears at low luminosities, 
but there are ample evidences that at least the broad line region (BLR) disappears at low 
luminosities. This happens for two known class of objects: 1) BL Lac sources and 2) low- 
luminosity type-2 AGN. 

BL Lacs are a class of (low-luminosity) blazars characterized by the absence of emission 
lines in their optical spectrum. Their broad-ban d spectrum is normal ly well understood in 
terms of the Synchrotron- Self Compton model (IMaraschi et al.lll992l ) where the electrons 
responsible for the synchrotron emission are up-scattering (via inverse Compton) the same 
synchrotron radiation to high energy. The main difference with the more luminous flat spec- 
trum radio quasars is that in these latter ones an additional high-energy component (refereed 
to as 'external Compton' component) is normally detected. In the external Compton model 
photons from the BLR and/or the disk are up-scattered to high energy by the electrons in 
the jet. It is believed that the absence of lines and of this external components is caused, in 
BL Lac objects, by the lack of the BLR, and a lower radiation field density. Moreover, BL 
Lacs are characterized by low Eddington ratios (e.g. \Edd — 0.01). Around this Eddington 
ratio, the accretion process experiences a transition from an optically thic k, geometrically 
thin disk (A^;^^ > 0.01) to a radiatively inefficient geometrically thick disk (iGhisellini et al. 
2OO9I ). Therefore in blazars with low Eddington ratios, the electrons in the jet find a medium 
starved of external radiation, weak or no lines are produced, and the AGN is classified as 
BL Lac. 

Another evidence for the absence of BLR at low luminosities is produced by low- 
luminosity type-2 AGN. These AGN, are sometimes referred to a s 'true' ty pe-2 AGN because 
when observed in polarized light they do not show broad lines (IHol 120081 ) and hence lack a 
BLR. How the "local" changes in the accretion regime aff ect the environm ent at the BLR and 



the to rus region is non-trivial and debated (according to iNicastrd (120001 ) and iNicastro et al. 



( I2OO3I ). the main driver for the disappearance of the BLR appears to be the accretion rate), 
nonetheless in this framework we find a convincing interpretation of our findings. The 
clumpy torus, and the BLR progressively disappear, resulting in a less efficient X-ray ob- 
scuration. Therefore we should expect F% to flatten or even to invert its dependence on the 
l uminosity at progres sively low Edd ington ratios. With f uture X-ray missions like NuSTAR 
( iHarrison et al.l 120051 ) . and NHXM (jPareschi et al.ll2009l ) it will be possible to sample with 
better statistics the population of the low-luminosity AGN and to investigate the behavior 
of the BLR and the torus in greater detail. 
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8. Summary and conclusion 

We addressed the study of a complete, flux limited, sample of local AGN collected by 
the Swift-BAT instrument in the flrst three years of survey. The sources are listed in Tab.O 
along with their properties. The aim of this work is to characterize the AGN population from 
two fundamental observables such as the hard X-ray (15-55 keV) luminosity and absorption. 
To this aim we jointly fltted the BAT spectra with the available follow up in the 0.3-10 keV 
domain. In the following we briefly review the main flndings. We remind the reader that 
AGN are deflned 'absorbed' if the column density for photoelectric absorption exceeds 10^^ 
atoms cm~^. 

• Performing a stacked analysis of the complete 199 AGN sample, a simple power law 
model was shown not to account for the continuum emission. In addition we performed 
the stacked analysis of the different subsamples of sources: unabsorbed, absorbed, and 
-for the flrst time- Compton-thick one. The average spectrum of CT sources was found 
to be dominated (in the BAT band) by the scattered component and its photon index 
was found to be compatible with ~2.0. According to our results, only ~30% of the 
source intrinsic flux is visible in the 15-55 keV band. 

• We showed that absorbed AGN are characterized by slightly harder spectra (1.91) with 
respect to the unabsorbed ones (2.00). Nonetheless the distributions are quite broad, 
resulting in a Kolmogorov-Smirnov probability of 3.5 x 10"^ of belonging to the same 
parent population. 

• We computed the observed Nh distribution, which shows that the observed fraction 
absorbed sources is 54±4%, with columns peaking at 10^^ cm^^. The observed fraction 
of Compton-thick objects is 4.6^^5%, a factor 2 to 3 lower than what predicted from 
population synthesis models at the fluxes of the BAT survey. 

• We estimated the bias of the BAT instrument against the detection of the Compton- 
thick AGN. We consequently derived the intrinsic Nh distribution by integrating the 
logN-logS in bins of logNn setting the minimum observed flux of integration so that 
the limiting intrinsic flux was the same for all the bins. Therefore we showed that even 
if the CT objects are only a minor fraction of the observed sample, their contribution 
rises to 20lg % in the intrinsic AGN population. 

• The relation between the observed fraction of obscured AGN and the hard X-ray lumi- 
nosity (-FVo), was found to have different behaviors according to the luminosity regime 
considered. For luminosities greater then 10^^ erg s~^ we found a monotonic decline 
with a slope of -0.26±0.05. At smaller luminosities, albeit affected by poor statistics. 
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we found a flattening of F% which we interpreted as the manifest disappearance of the 
obscuring region. In a disk-cloud outflow scenario, this is indeed expected to happen 
under a critical luminosity, which is of the same order of the luminosity at which we 
observe the flattening. 

• We showed that the obscuration-luminosity relation can be explained by the different 
X-ray Luminosity Functions of the obscured and unobscured subsamples. This in turn 
means that absorbed AGN are intrinsically less luminous. This result, if the mass 
distribution is narrow, points towards a trend in the Eddington ratios in which objects 
accreting at lower values (thus having smaller effects on their environments) are more 
absorbed. Obscuring clouds would be able to come closer to the nuclear region without 
being affected, and bury the AGN. 



A key test to improve our findings will be the calculation of the BH masses for the 
sample. By means of the two physical quantities (i.e. A^;^^ and mass) we plan to test 
whether it is possible to find a sequence relating absorption, Eddington ratio , and mass 
of the black hole . Sorn e information can already be found in iMiddleton et al.l ( 120081 ) and 



Beckmann et al.l (120091 ). They found that for two different samples of AGN selected in the 



local Universe, the mean Eddington ratio is in the 0.01-0.06 range. ICappelluti et al.l (|2010[ ) 
also showed tha t, on average, the BAT AG N used in this work have an Eddington ratio of 
0.01. Moreover, IMiddleton et al.l (120081 ) and iBeckmann et al.l (120091 ) found that on average 
unobscured AGN have larger Eddington ratios with respect to obscured ones. This would 
be consistent with the presence of a trend in \Edd for non-jetted AGN. The obscured objects 
could be accreting at lower Eddington ratios and with flatter spectra; the unobscured ones, 
with steeper spectra, could be accreting slightly more efficiently. By means of this analysis it 
would be also possible to relate our findings to a physical consistent picture, in the framework 
of merger driven AGN activity. Indeed, merging of gas-rich galaxies provides an efficient way 
to fu nnel large amount of gas and dust to the central black hole and trigger s AGN activity 



(e.g. 



Kauffmann fc Haehneltl 120001: IWvithe fc Loebll2003l : ICroton et al.ll2006t l. Noteworthy a 
very recent work by iKoss et al.l ( 120101 ) determined that a considerable fraction (i.e. ~24 %) 
of the hosts of the BAT AGN have a close companion within 30 kpc and are experiencing a 
major merging event. This fraction is extremely relevant when compared to a control sample 
of local (i.e. 2; < 0.1) optically selected narrow-line AGN, where the fraction of interacting 
companion s is ~ 1%. We counted how many objects of the sample of merging AGN of 
Koss et al.l ( 120101 ) are obscured. We found that >63 % are obscured AGn| with an average 
column density of logNn ~ 23.4. A limitation to this simple picture is that in case of a 



^ For comparison, we remind that in our sample the fraction of obscured sources is ^50 %. 
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merging event and of a large gas quantity being funne led towards the center, the black hole 
is exp ected to accre te with high E ddington ratios (e.g. iDi Matteo et al.ll2005l : iHopkins et al. 



20061 ). According to iFabiaru ( 119991 ) a SMBH accreting at Eddington luminosities should clear 
the environment from any Compton-thin (e.g. logNj:/ < 24) column density and therefore 
transit to a less obscured phase. If the very first phase of the gas-rich merger event is the 
creation of a Compton-thick AGN, then we would expect it to display large Eddington ratios. 
However, two of the most famous Compton-thick AGN (Circinus and N GC 4945) display 



an Eddington ratio far from unity (i.e. < 10 ^, see iGiiltekin et al.l 120091 ). proving that this 



argument still escapes a conclusive explanation. The largest hard X-ray selected samples of 
AGN in the local Universe may shed some light on the physical interpretation of the feedback 
of black holes on their surroundings. 
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Table 5. Swift sample of AGN. 
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23.6 


(1) 


J0342. 2-2114 


55.5.54 


-21.244 


2.15±0.18 


11.8 


SWIFT J0342.0-2115 


Syl 


0.01 


-1 OO+0.20 


20.5 


(1) 


J0347.3-3029 


56.850 


-30.500 


0.89±0.17 


5.3 


RBS 0741 


Syl 


0.10 


2.04+111 


20.4 


(2) 


J0350.7-5022 


57.679 


-50.377 


1.29±0.17 


7.5 


SWIFT J0350. 1-5019 


Sy2 


0.04 


1.79+°-3| 


23.2 


(2) 


J0357.0-4039 


59.268 


-40.666 


0.89±0.17 


5.4 


2MASX J03565655-4041453 


Syl.9 


0.07 


1.56tJ3-l 


22.5 


(1) 


J0402.5-1804 


60.639 


-18.077 


1.35±0.19 


7.0 


ESO 549- G049 


Sy2 


0.03 


-, 7C+0.50 
38 


22.4 


(2) 


J0407. 5+0342 


61.883 


3.717 


1.90±0.25 


7.6 


3C 105 


Sy2 


0.09 


1-91+^28 


23.4 


(1) 


J0415.2-0753 


63.800 


-7.900 


1.31±0.23 


5.6 


LEDA 14727 


Syl 


0.04 


9 9t--f0.44 

^•^O Q 


23.5 


(2) 


J0426.4-5712 


66.603 


-57.201 


1.40±0.17 


8.2 


IH 0419-577 


Syl 


0.10 


9 t;7+0:66 
^•o 1 Q 22 


19.5 


(1) 


J0433. 4+0521 


68.355 


5.365 


5.21±0.26 


19.8 


3C-120 


Syl 


0.03 


9 19-1-0:13 

12 


21.2 


(1) 


J0438. 5-1049 


69.633 


-10.830 


1.48±0.23 


6.4 


MCG-02-12-050 


Syl 


0.04 


2.02t»:6i 


20.0 


(2) 


J0444. 7-2812 


71.199 


-28.200 


1.07±0.18 


5.9 


2MASX J04450628-2820284 


Sy2 


0.15 


2 x7+"'^^ 


20.0 


(2) 


J0451. 8-5807 


72.966 


-58.133 


0.88±0.17 


5.2 


RBS 0594 


Syl 


0.09 




20.0 


(2) 


J0453.5+0403 


73.380 


4.060 


2.11±0.28 


7.6 


CGCG 420-015 


Sy2 


0.03 


2.04l°i 


~24.2 


(4) 


J0455. 3-7528 


73.841 


-75.477 


1.27+0.18 


6.9 


ESO 33-2 


Sy2 


0.02 


9 C9-I-O.49 
''—0.41 


22.1 


(2) 


J0505. 9-2351 


76.497 


-23.854 


2.78+0.20 


13.9 


XSS J05054-2348 


Sy2 


0.04 




22.7 


(1) 


J0516. 2-0009 


79.071 


-0.161 


4.11+0.28 


14.5 


QSO B0513-002 


Syl 


0.03 




20.0 


(2) 


J0519.7-3240 


79.930 


-32.676 


2.38+0.19 


12.9 


SWIFT J0519.5-3140 


Sy2 


0.04 


1 79-1-0.18 
-^•'^-0.14 


21.1 


(1) 


J0519.8-4546 


79.963 


-45.774 


2.49±0.17 


14.9 


Pictor-A 


Syl 


0.04 




21.0 


(1) 


J0524.2-1212 


81.050 


-12.200 


1.42+0.25 


5.6 


LEDA 17233 


Syl 


0.05 




21.2 


(2) 


J0552. 3-0727 


88.090 


-7.457 


14.75+0.29 


51.7 


NGC 2110 


Sy2 


0.01 


'-■'^-0.00 


22.6 


(1) 


J0552.3+5929 


88.100 


59.500 


1.14+0.21 


5.3 


IRAS 05480+5927 


Syl 


0.06 


44+1.00 


21.1 


(2) 


J0558. 1-3820 


89.549 


-38.347 


2.12±0.18 


11.6 


EXO 055620-3820.2 


Syl 


0.03 


9 9i-f0.23 
z.zi_Q 21 


22.2 


(1) 


J0602.9-8633 


90.749 


-86.555 


1.82±0.22 


8.4 


SWIFT J0601.9-8636 


Sy2 


0.01 


, ^,7-1-0.24 
^■"'-0.37 


~24 


(4) 


J0603. 1+6523 


90.799 


65.399 


1.38+0.20 


6.8 


UGC 3386 


GALAXY 


0.02 


^■J^'J-0.33 


23.2 


(3) 


J0615. 8+7101 


93.967 


71.021 


6.08+0.20 


30.8 


Mrk 3 


Sy2 


0.01 




24.0 


(1) 


J0623. 9-3214 


95.994 


-32.248 


1.53+0.20 


7.5 


ESO 426-G 002 


Sy2 


0.02 




23.9 


(2) 


J0624. 1-6059 


96.028 


-60.998 


1.25±0.17 


7.4 


SWIFT J2141.0+1603 


Sy2 


0.04 


9 CI +0.42 


23.4 


(2) 


J0640.7-4324 


100.200 


-43.400 


0.92±0.18 


5.2 


2MASX J06400609-4327591 


Sy2 


0.06 


Z.UD_Q gg 


23.4 


(2) 
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SWIFT NAME R.A. 


DEC 


Flux [15-55 keV] S/N 


ID 


Type* 


redshift Ph. Index log(NH) Reference 


(J2000) 


(J2000) 


(10-11 Cgs) 






(BAT only) 



J0652. 1+7425 


103.044 


74.425 


3.29+0.19 


17.1 


Mrk 6 


Syl.5 


0.02 




23.0 


(1) 


J0656.1+3959 


104.027 


39.986 


2.29+0.26 


8.7 


UGC 3601 


Syl 


0.02 


1.97l°i9 


21.3 


(2) 


J0718.0+4405 


109.517 


44.084 


1.67+0.24 


7.1 


2MASX J07180060+4405271 


Syl 


0.06 


r, 9,-1-0.41 


20.0 


(2) 


J0742.5+4947 


115.644 


49.793 


2.98+0.21 


14.4 


Mrk 79 


Syl.2 


0.02 


— u. lb 


20.8 


(1) 


J0800.1+2322 


120.032 


23.370 


1.62+0.24 


6.6 


SDSS J0759.87+232448.3 


GALAXY 


0.03 


1.70+°|J 


22.3 


(2) 


J0800.3+2638 


120.099 


26.648 


1.79+0.24 


7.5 


IC 486 


Syl 


0.03 


1.8010^0 


22.2 


(2) 


J0804.2+0506 


121.050 


5.101 


3.18+0.25 


12.6 


UGC 4203 


Sy2 


0.01 


9 t;o+0.60 


23.5 


(10) 


J0811. 1+7602 


122.798 


76.049 


1.15+0.19 


6.2 


PC 0804+761 


Syl 


0.10 


2.58+°i° 


20.0 


(2) 


J0814.4+0423 


123.600 


4.400 


1.26+0.24 


5.2 


CGCG 031-072 


Syl 


0.03 


i-ao+oS 


23.3 


(2) 


J0823.2-0456 


125.800 


-4.947 


1.35+0.23 


6.0 


SWIFT J0823.4-0457 


Sy2 


0.02 


1.66+°-38 


23.5 


(3) 


J0832. 8+3706 


128.200 


37.100 


1.03+0.20 


5.2 


RBS 707 


Syl.2 


0.09 


2.41+«/g 


20.0 


(2) 


J0839.8-1214 


129.950 


-12.248 


1.28+0.21 


6.1 


3C 206 


Syl 


0.20 


2.04+°:|t 


21.0 


(2) 


J0904.9+5537 


136.250 


55.632 


1.04+0.17 


6.0 


SWIFT J0904.3+5538 


Syl.5 


0.04 




21.0 


(1) 


J0911.5+4528 


137.898 


45.471 


1.23+0.18 


7.0 


SWIFT J0911.2+4533 


Sy2 


0.03 


2.45+°it 


23.5 


(1) 


J0918.4+1618 


139.615 


16.316 


1.65+0.21 


8.0 


Mrk 704 


Syl.5 


0.03 


1.98l«-26 


21.5 


(2) 


J0921.0-0803 


140.257 


-8.067 


2.59+0.20 


12.9 


SWIFT J0920.8-0805 


Sy2 


0.02 


9 1 c; + 0.21 


22.8 


(1) 


J0923.8+2256 


140.962 


22.936 


2.05+0.20 


10.5 


MCG +04-22-042 


Syl.2 


0.03 


i-oO Q 21 


20.6 


(1) 


J0925.2+5217 


141.316 


52.285 


3.03+0.17 


18.1 


Mrk 110 


Syl 


0.04 




20.6 


(1) 


J0945.8-1419 


146.468 


-14.332 


1.28+0.21 


6.1 


NGC 2992 


Sy2 


0.01 


1>55 Q 22 


22.0 


(1) 


J0947.7-3056 


146.939 


-30.948 


11.50+0.22 


51.5 


ESQ 434-40 


Sy2 


0.01 


9 r,7+o:oi 

^•^ ' _() 01 


22.2 


(3) 


J0959.6-2250 


149.916 


-22.834 


4.44+0.22 


19.8 


NGC 3081 


Sy2 


0.01 


i.8ol«:i« 


23.5 


(1) 


JlOOl. 8+5542 


150.453 


55.700 


1.43+0.16 


8.8 


NGC 3079 


Sy2 


0.004 


-I OO+0.26 
^•*°-0.25 


24.7 


(4) 


J1006.0-2306 


151.500 


-23.100 


1.25+0.23 


5.5 


ESQ 499-G 041 


Syl 


0.01 


1 p-Y+0.60 
-0.62 


21.4 


(2) 


J1021. 7-0327 


155.450 


-3.450 


1.35+0.22 


6.3 


MCG+00-27-002 


Syl 


0.04 


9 04+0.54 
^••5^-0.44 


20.0 


(2) 


J1023.5+1951 


155.888 


19.864 


7.35+0.20 


36.8 


NGC 3227 


Syl.5 


0.004 




22.8 


(1) 


J1031.8-3451 


157.975 


-34.860 


4.71+0.26 


18.4 


NGC 3281 


Sy2 


0.01 




23.9 


(6) 


J1031.9-1417 


157.996 


-14.300 


2.13+0.23 


9.3 


H 1029-140 


Syl 


0.09 


9 -,7+0.32 
-0.29 


20.0 


(2) 


J1044.0+7023 


161.003 


70.400 


0.99+0.17 


5.9 


MCG+12-10-067 


Sy2 


0.03 


i.DU_Q 35 


23.3 


(2) 


J1046.5+2556 


161.649 


25.950 


1.12+0.19 


5.9 


UGC 05881 


GALAXY 


0.02 


1 60+°-^* 

i.DU_g 3g 


23.0 


(2) 


J1048.5-2512 


162.149 


-25.200 


1.42+0.27 


5.3 


NGC 3393 


Sy2 


0.01 


9 (,,+0.43 


24.7 


(4) 


J1049.3+2256 


162.350 


22.950 


1.57+0.20 


8.0 


SWIFT J1049.4+2258 


Sy2 


0.03 


, og+0.24 

^•S9-0.23 


23.3 


(2) 


J1106.6+7234 


166.654 


72.571 


6.45+0.17 


38.0 


NGC 3516 


Syl.5 


0.01 


1.901H1 


21.2 


(1) 


J1115.9+5426 


168.999 


54.450 


0.88+0.15 


5.7 


SDSS J111519.98+542316.6 


Sy2 


0.07 


1.80l°;S 


22.4 


(2) 


J1125.4+5421 


171.352 


54.351 


0.97+0.15 


6.3 


ARP 151 


Syl 


0.02 


-, 79+0.48 


20.0 


(2) 
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SWIFT NAME R.A. 


DEC 


Flux [15-55 keV] S/N 


ID 


Type* 


redshift Ph. Index log(NH) Reference 


(J2000) 


(J2000) 


(10-11 Cgs) 






(BAT only) 



J1127.5+1908 


171.900 


19.148 


1.12+0.20 


5.6 


IRXS J112716.6+ 190914 


Syl 


0.10 


1.84l«-g 


21.4 


(2) 


J1132. 7+5259 


173.188 


52.988 


1.01+0.15 


6.6 


UGC 6527 


Syl 


0.03 


1 oc: + 0.36 
i.oO Q 34 


20.0 


(3) 


J1136.5+2132 


174.150 


21.548 


1.13+0.19 


5.9 


Mrk 739 


Syl 


0.03 




20.7 


(2) 


J1139.0-3744 


174.764 


-37.741 


10.07+0.27 


37.9 


NGC 3783 


Syl 


0.01 


1.94+0-01 


22.5 


(1) 


J1139. 1+5912 


174.783 


59.212 


1.25+0.15 


8.1 


SBS 1136+594 


Syl. 5 


0.06 




19.6 


(1) 


J1139.4+3156 


174.869 


31.935 


1.00+0.17 


5.8 


NGC 3786 


Syl. 8 


0.01 


1.74+«-« 


22.5 


(3) 


J1144.7+7939 


176.190 


79.662 


2.13+0.18 


11.9 


SWIFT J1143.7+7942 


Syl. 2 


0.02 




20.6 


(1) 


J1145.3+5859 


176.349 


59.000 


0.81+0.15 


5.3 


Ark 320 


GALAXY 


0.01 


9 90+0.74 

^•^O Q 


22.2 


(2) 


J1145.5-1825 


176.393 


-18.428 


2.84+0.27 


10.5 


2MASX J 11454045- 1827149 


Syl 


0.03 


2.02+^-2? 


20.5 


(1) 


J1148. 9+2938 


177.230 


29.634 


1.06+0.18 


6.1 


MCG+05-28-032 


LINER 


0.02 


1.85+»-3,?, 


22.5 


(2) 


J1158. 0+5526 


179.502 


55.449 


1.04+0.15 


6.9 


NGC 3998 


Syl 


0.003 


2.05+i]l'; 


20.1 


(8) 


J1201.0+0647 


180.250 


6.800 


1.18+0.21 


5.6 


SWIFT J1200.8+0650 


Sy2 


0.04 


1.89+«-i 


22.8 


(1) 


J1203.0+4432 


180.773 


44.534 


2.33+0.15 


15.1 


NGC 4051 


Syl.5 


0.002 


2.45+°-;f 


20.5 


(1) 


J1204.5+2018 


181.149 


20.301 


1.32+0.19 


7.1 


ARK 347 


Sy2 


0.02 


1.76+O-i 


23.2 


(1) 


J1206. 2+5242 


181.565 


52.710 


1.24+0.15 


8.3 


NGC 4102 


LINER 


0.003 


1.74+0-30 


20.9 


(1) 


J1209. 1+4700 


182.300 


47.000 


0.76+0.15 


5.0 


Mrk 198 


Sy2 


0.02 


1-76+21^ 


22.8 


(2) 


J1209. 4+4341 


182.370 


43.686 


1.56+0.15 


10.1 


NGC 4138 


Syl.9 


0.003 




22.9 


(1) 


J1210.5+3924 


182.633 


39.406 


24.60+0.16 


153.4 


NGC 4151 


Syl.5 


0.003 


1 QQ+0.01 

'-^'-'—0 01 


22.5 


(1) 


J1210.6+3819 


182.667 


38.333 


0.96+0.16 


6.0 


LEDA 38759 


Syl 


0.02 


, ni -1-0:36 
— 33 


22.6 


(2) 


J1217.2+0711 


184.300 


7.200 


1.21+0.20 


5.9 


NGC 4235 


Syl 


0.01 


1 Kq+0:44 


21.2 


(3) 


J1218. 3+2950 


184.593 


29.839 


1.53+0.17 


9.1 


Mrk 766 


Syl.5 


0.01 




21.7 


(1) 


J1219.0+4715 


184.750 


47.252 


0.96+0.15 


6.3 


NGC 4258 


Syl 


0.001 


J^-^^-0.34 


22.9 


(11) 


J1222.0+7518 


185.503 


75.311 


1.27+0.17 


7.4 


Mrk 205 


Syl 


0.07 


9 co+0.42 


20.7 


(12) 


J1225.7+1239 


186.447 


12.665 


12.58+0.19 


65.6 


NGC 4388 


Sy2 


0.01 




23.6 


(1) 


J1225. 8+3330 


186.466 


33.513 


1.25+0.16 


7.7 


NGC 4395 


Syl 


0.001 


9 1 r; + 0.27 


22.3 


(1) 


J1235. 6-3955 


188.902 


-39.919 


10.21+0.26 


39.4 


NGC 4507 


Sy2 


0.01 


1.981";- 


23.5 


(1) 


J1238. 8-2718 


189.723 


-27.308 


4.39+0.28 


15.9 


ESQ 506-027 


Sy2 


0.02 


1.74±«-3 


23.6 


(9) 


J1239.0-1611 


189.769 


-16.196 


2.02+0.26 


7.7 


XSS J12389-1614 


Sy2 


0.04 


l.b»_o 22 


22.5 


(1) 


J1239.5-0520 


189.898 


-5.341 


4.52+0.23 


19.8 


NGC 4593 


Syl 


0.01 


1 Q1-l-0.il 

J^-yj^-0.11 


20.3 


(1) 


J1246. 6+5434 


191.661 


54.575 


1.34+0.15 


9.0 


NGC 4686 


LINER 


0.02 


, 7t-+0.25 


23.8 


(2) 


J1302. 8+1624 


195.700 


16.400 


0.90+0.17 


5.1 


Mrk 0783 


Syl. 2 


0.07 


1 Q2+"'*^ 
-'-■■^^-().41 


21.0 


(2) 


J1306. 7-4024 


196.698 


-40.415 


2.37+0.27 


8.9 


ESQ 323-077 


Syl. 2 


0.02 


2 n3+"'^^ 


22.7 


(2) 


J1309.1+1137 


197.279 


11.632 


2.19+0.18 


12.0 


SWIFT J1309.2+1139 


Sy2 


0.03 


l-Do_Q J5 


23.4 


(13) 


J1315.4+4424 


198.852 


44.404 


1.28+0.15 


8.4 


IGR J13149+4422 


Sy 


0.04 


9 90-1-0.31 
^•^o_0.28 


22.8 


(2) 



CO 
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SWIFT NAME R.A. 


DEC 


Flux [15-55 keV] S/N 


ID 


Type* 


redshift Ph. Index log(NH) Reference 


(J2000) 


(J2000) 


(10-11 Cgs) 






(BAT only) 



J1322.3-1642 


200.591 


-16.716 


2.57±0.27 


9.5 


MCG -03-34-064 


Syl.8 


0.02 


9 1C+0.30 

—0 28 


23.6 


(1) 


J1325.4-4301 


201.366 


-43.017 


49.77±0.26 


187.6 


Cen A 


Sy2 


0.002 


1 ot; + 0:00 

—0 00 


22.7 


(1) 


J1334.8-2323 


203.700 


-23.400 


1.49±0.29 


5.1 


ESQ 509-38 


Sy2 


0.03 


54 


20.0 


(2) 


J1335.7-3418 


203.944 


-34.302 


4.86±0.29 


16.6 


MCG -06-30-015 


Syl.2 


0.01 


2.24+°-';? 


21.7 


(1) 


J1338.1+0133 


204.547 


4.552 


3.63±0.20 


17.9 


NGC 5252 


Sy2 


0.02 


1.671S-12 


22.64 


(13) 


J1341.4+3022 


205.356 


30.369 


1.15±0.16 


7.2 


Mrk 268 


Sy2 


0.04 


9 OO + 0.35 


23.3 


(2) 


J1349.5-3018 


207.390 


-30.304 


17.87±0.31 


58.2 


IC 4329A 


Syl 


0.02 


2.05+»:»» 


21.6 


(1) 


J1353.2-I-6919 


208.305 


69.327 


2.78±0.17 


16.7 


Mrk 279 


Syl.5 


0.03 




20.5 


(1) 


J1356.1-I-3835 


209.033 


38.583 


1.22±0.16 


7.7 


Mrk 464 


Syl 


0.05 


1.69+°-30 


20.0 


(3) 


J1408.4-3024 


212.100 


-30.400 


1.65±0.32 


5.1 


PGC 050427 


Syl 


0.02 


2.67lJ5-«^ 


21.2 


(2) 


J1413.5-0312 


213.375 


-3.201 


14.38±0.24 


59.0 


NGC 5506 


Syl.9 


0.01 


2 27+('-ii 


22.5 


(1) 


J1418.2+2507 


214.568 


25.133 


3.12±0.17 


18.2 


NGC 5548 


Syl.5 


0.02 




20.4 


(1) 


J1419.5-2639 


214.893 


-26.663 


3.49±0.34 


10.4 


ESQ 511-G030 


Syl 


0.02 


9 11+0.22 

J- -L Q 2Q 


21.2 


(1) 


J1421.6-I-4750 


215.420 


47.838 


1.03±0.16 


6.4 


QSO B1419-I-480 


Syl 


0.07 


1 7o-H0:36 
J-' 1 <-> Q 27 


21.3 


(3) 


J1424.3+2435 


216.100 


24.600 


0.89±0.17 


5.1 


NGC 5610 


GALAXY 


0.02 


1 

i.oo Q 42 


22.8 


(2) 


J1429.6-I-0117 


217.400 


1.300 


1.26±0.23 


5.4 


QSO B1426-I-015 


Syl 


0.09 




20.0 


(3) 


J1436.5-I-5847 


219.149 


58.798 


1.37±0.16 


8.3 


QSO J1436-f5847 


Syl 


0.03 




23.5 


(2) 


J1441.2-I-5330 


220.300 


53.500 


0.85±0.16 


5.1 


Mrk 477 


Sy2 


0.04 


1.47+°-g 


24.0 


(14) 


J1442.6-1713 


220.664 


-17.223 


4.83±0.34 


14.4 


NGC 5728 


Sy2 


0.01 


i.84+o-';t 


24.3 


(4) 


J1453.1+2556 


223.282 


25.936 


1.29±0.18 


7.0 


RX J1453.1-I-2554 


Syl 


0.05 


1.8210^6 


20.0 


(2) 


J1504.2+1025 


226.073 


10.417 


1.51±0.22 


6.8 


Mrk 841 


Syl 


0.04 




21.3 


(1) 


J1515.4-I-4201 


228.868 


42.033 


1.05±0.18 


5.9 


NGC 5899 


Sy2 


0.01 


-1 79 + 0.56 
-^•'^-0.51 


23.1 


(2) 


J1536.2-I-5753 


234.061 


57.890 


1.43±0.18 


7.9 


Mrk 290 


Syl 


0.03 


9 1Q+0.39 
^•J^y_0.34 


20.4 


(1) 


J1548.4-1344 


237.106 


-13.749 


2.91±0.39 


7.4 


NGC 5995 


Sy2 


0.03 


2 01+0-26 
^•<-'i_o.24 


22.0 


(2) 


J1554.8+3242 


238.700 


32.700 


1.04±0.20 


5.2 


2MASX J15541741+3238381 


Syl 


0.05 


1 66+"'^*' 


20.0 


(2) 


J1628.3+5147 


247.082 


51.793 


2.45±0.20 


12.4 


SWIFT J1628.1-t-5145 


Syl.9 


0.05 


9 1 O + 0.19 


23.3 


(1) 


J1653.2+0224 


253.319 


2.404 


4.20d=0.35 


12.1 


NGC 6240 


Sy2 


0.02 


2.09t«:i« 


24.3 


(4) 


J1822.1-I-6421 


275.541 


64.361 


1.10±0.21 


5.3 


QSO B1821-I-643 


Syl 


0.30 


2.9lir6^ 


20.0 


(15) 


J1824.2-5620 


276.057 


-56.348 


1.98±0.29 


6.9 


IC 4709 


Sy2 


0.02 


1 Q4+0-28 


23.1 


(2) 


J1835. 1+3240 


278.791 


32.683 


4.67±0.21 


22.0 


3C 382 


Syl 


0.06 


2 nq+0.11 


21.1 


(1) 


J1837.1-5922 


279.284 


-59.368 


1.79±0.28 


6.3 


FAIRALL 49 


Sy2 


0.02 


9 t:t; + 0.43 


22.3 


(14) 


J1838.6-6523 


279.658 


-65.394 


6.23±0.28 


22.4 


ESQ 103-035 


Sy2 


0.01 


9 1S+0-12 
^•-'-*-0.12 


23.2 


(1) 


J1842.4+7946 


280.616 


79.771 


5.82±0.19 


29.9 


3C 390.3 


Syl 


0.06 


1 Qp; + 0.07 

l.»0_Q 


21.0 


(1) 


J1845. 1-6223 


281.297 


-62.399 


2.52±0.28 


9.0 


ESQ 140-43 


Syl 


0.01 


9 17+0.27 
^■^'-0.25 


22.4 


(2) 



Table 5 — Continued 



SWIFT NAME R.A. 


DEC 


Flux [15-55 kcV] S/N 


ID 


Type"!" redshift Ph. Index log(NH) Reference 


(J2000) 


(J2000) 


(10-11 cgs) 




(BAT only) 



J1857.3-7827 


284.341 


-78.464 


1.86±0.26 


7.1 


LEDA 140831 


Syl 


0.04 


1 


otr-l-0.36 

— 34 


20.0 


(2) 


J1921. 2-5840 


290.323 


-58.677 


3.54±0.28 


12.6 


ESQ 141-55 


Syl 


0.04 


1 


00+0.19 

18 


20.0 


(3) 


J1942. 7-1018 


295.680 


-10.316 


4.30±0.30 


14.3 


NGC 6814 


Syl 


0.01 


1 


qi-l-O.lS 

— 14 


20.8 


(1) 


J2009.1-6103 


302.289 


-61.064 


3.03±0.26 


11.5 


SWIFT J2009.0-6103 


Syl 


0.01 


1 




21.8 


(1) 


J2044. 1-1043 


311.039 


-10.731 


5.62±0.29 


19.7 


Mrk 509 


Syl. 2 


0.03 


2 


04+0-}? 


20.7 


(1) 


J2052.0-5703 


313.017 


-57.063 


4.63±0.25 


18.3 


IC 5063 


Sy2 


0.01 


1 


89tO-}^ 


23.3 


(1) 


J2109.1-0939 


317.300 


-9.652 


1.58±0.27 


5.9 


IH 2107-097 


Syl 


0.03 


2 


1 t:+0.40 
— 35 






J2132. 1-3343 


323.028 


-33.727 


2.92±0.27 


10.7 


CTS 109 


Syl 


0.03 


2 


07+0.22 

— 20 


20.0 


(2) 


J2136.0-6223 


324.006 


-62.400 


2.42±0.23 


10.6 


QSO J2136-6224 


Syl 


0.06 


2 


Qo-f0.'26 

24 


20.0 


(2) 


J2138.8+3206 


324.713 


32.115 


1.29±0.00 


6.3 


LEDA 67084 


Syl 


0.03 


2 




20.0 


(2) 


J2200.7-fl033 


330.199 


10.565 


1.76±0.21 


8.6 


SWIFT J2200.9-I-1032 


Syl. 9 


0.03 


2 




22.2 


(1) 


J2202. 1-3152 


330.526 


-31.878 


8.01±0.25 


31.8 


NGC 7172 


Sy2 


0.01 


1 




22.9 


(1) 


J2204.5-I-0335 


331.149 


3.600 


1.33±0.21 


6.3 


IRAS 22017-1-0319 


Sy2 


0.06 


2 


nQ+0.48 

40 


22.5 


(2) 


J2209.5-4709 


332.387 


-47.166 


3.02±0.23 


13.4 


NGC 7213 


Syl. 5 


0.03 


1 


q2+0.l7 

16 


20.6 


(1) 


J2223.8-0207 


335.962 


-2.121 


1.93±0.22 


8.9 


3C 445 


Syl 


0.06 


1 


qq+0-'24 


23.2 


(2) 


J2235.8-2603 


338.966 


-26.054 


2.76±0.24 


11.6 


NGC 7314 


Syl.9 


0.005 


2 


04+0-20 


21.8 


(1) 


J2236. 1-1-3357 


339.040 


33.952 


1.66±0.19 


8.8 


Arp 319 


Sy2 


0.02 


1 


7t:+0.27 
'3-0.23 


23.7 


(2) 


J2236.8-1235 


339.223 


-12.599 


1.43±0.23 


6.2 


Mrk 915 


Syl 


0.02 


1 


c:q+0.33 
"J ="-0.27 


22.8 


(2) 


J2245. 7-1-3941 


341.449 


39.695 


1.71±0.18 


9.2 


3C 452 


Sy2 


0.08 


1 




23.4 


(1) 


J2254.1-1734 


343.535 


-17.578 


5.67±0.23 


24.7 


MR 2251-178 


Syl 


0.06 


2 


06+0-" 


20.8 


(1) 


J2258.9-I-4053 


344.749 


40.899 


1.31±0.18 


7.2 


UGC 12282 


Syl 


0.02 


1 


fit:+0.28 
"■^-0.22 


23.9 


(2) 


J2259.5-I-2455 


344.899 


24.929 


1.51±0.19 


8.0 


LEDA 70195 


Syl 


0.03 


1 


^"-0.28 


20.0 


(2) 


J2303.2-I-0853 


345.809 


8.885 


3.87±0.20 


19.4 


NGC 7469 


Syl.2 


0.02 


2 


14+0-12 
1^-0.12 


20.6 


(1) 


J2304.7-0841 


346.194 


-8.686 


6.09±0.22 


27.9 


Mrk 926 


Syl.5 


0.05 


1 


q7+0.01 
y -0.01 


21.1 


(1) 


J2304.7-I-1217 


346.200 


12.300 


1.11±0.20 


5.6 


NGC 7479 


Sy2 


0.01 


1 


Qo+0.40 
y°-0.37 


23.6 


(3) 


J2318.4-4221 


349.614 


-42.360 


4.09±0.20 


20.7 


NGC 7582 


Sy2 


0.01 


1 


qc+0.22 
^■^-0.24 


24.1 


(4) 


J2319.0-I-0014 


349.762 


0.241 


2.82±0.00 


13.5 


NGC 7603 


Syl 


0.03 


2 


07+0.18 
'^'-0.17 


20.0 


(2) 


J2326.3-I-2154 


351.600 


21.900 


0.97±0.19 


5.1 


RBS 2005 


Syl 


0.12 


2 


nK+0.52 
'-''^-0.47 


20.3 


(2) 


J2342.0-I-3035 


355.500 


30.600 


1.08±0.19 


5.7 


UGC 12741 


Sy2 


0.02 


1 


R4+0-33 
°*-0.30 


23.7 


(2) 



Ol 



t The optical classification comes mainly from Tueller et al. (2008), Winter et al. (2009a'). IParisi et al.l l|2009h . lCusumano et all (|2O10h . SIMBAD, and NED. 
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iMatt et al.ll2009l: (lll lCappi et al.ll2006l : (12l |Page et al ll2005l : (13'l lWinter et al.ll2008l : (14'l lShu et al.ll2007l: (15) Ijimcnez-Bailon et al.ll2007l " 



